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The  research  described  in  this  report  was  prepared  by  Honeywell  Inc. , 
Minneapolis,  Minnesota  55413,  under  Air  Force  Contract  F33615-75-C-3046. 
It  was  initiated  under  the  AFFDL  Task  Number  82190221,  "Optimal  Control 
of  Flexible  Aircraft,  "  Project  Number  8219,  "  Stability  and  Control  of  Aero¬ 
space  Vehicles.  "  This  work  was  directed  by  the  Control  Criteria  Branch 
(FGC),  Flight  Control  Division  of  the  Air  Force  Flight  Dynamics  Laboratory 
and  was  administered  by  Mr.  Charles  R.  Stockdale  of  the  Control  Criteria 
Branch.  Special  thanks  to  Mr.  Robert  C.  Schwanz  of  FGC  and  Mr.  Gary 
Grimes  of  ASD/ADDP  for  their  continued  support  toward  this  contract. 

The  technical  work  reported  in  this  volume  was  conducted  by  the  Research 
Department  at  the  Systems  and  Research  Center  of  Honeywell  Inc.  Dr.  A. 

F.  Konar  was  the  Honeywell  Program  Manager  and  the  principal  investigator 
on  this  contract.  He  was  assisted  by  Mr.  C.  R.  Stone,  Dr.  J.  K.  Mahesh, 
and  Miss  M.  Hank.  This  report  covers  work  done  from  April  1975  to  April 
1976. 

The  work  under  this  contract  was  reported  in  three  volumes  entitled  "Active 
Control  Synthesis  for  Flexible  Vehicles." 
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Volume  II  KONPACT  Program  Listing 

Volume  III  KONPACT  Users  Manual 
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SECTION  I 


INTRODUCTION 


The  general  objective  of  this  study  is  to  develop  techniques  and  tools  necessary 
for  rapid  design  of  an  active  control  system  for  aircraft  with  lightly  damped 
structural  modes.  The  synthesis  techniques  provided  here  are  aimed  at 
reducing  the  engineering  man-hours  presently  required  for  a  flight  control 
system  design,  thus  effecting  a  cost  reduction.  Improvements  in  the  fatigue 
life,  ride  qualities,  and/or  handling  qualities  of  military  aircraft  are  sought 
by  controlling  the  lightly  damped  modes  and  thus  improving  their  mission 
performance. 

The  present  scope  of  this  study  is  to  develop  programs  to  interface  the  level 
2.01.00  FLEXSTAB  computer  program  system  with  existing  Air  Force- 
owned  optimal  control  computer  programs.  These  programs  represent  ad¬ 
vanced  computational  techniques  required  to  perform  quantitative  analysis 
of  multisurface  control  systems.  The  resulting  interface  program  system 
is  called  "KONPACT  -  Computer  Programs  for  Active  Control  Technology.  " 
KONPACT  provides  the  capability  to  model,  synthesize,  analyze,  and  design 
automatic  control  systems  efficiently  working  together  with  FLEXSTAB.  It 
can  also  be  used  as  a  stand  alone  program. 

The  work  performed  under  this  contract  is  reported  in  three  volumes: 

Volume  I.  KONPACT  Theoretical  Description  and  Demonstration 
Volume  II.  KONPACT  Program  Listing 
Volume  III.  KONPACT  Users  Manual 
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This  document  reports  the  analytical  techniques  and  algorithms  used  in 
KONPACT.  It  also  demonstrates  how  these  techniques  are  applied  to 
flexible  aircraft  control  system  design. 

An  overview  of  the  Honeywell  theoretical  approach  to  control  system  design 
using  KONPACT  is  provided  in  Section  II.  The  process  of  the  Active  Lift 
Distribution  Control  System  (ALDCS)  design,  brief  description  of  KONPACT, 
and  the  application  results  for  the  C-5A  vehicle  are  described  in  this  section. 

In  Section  III  the  technique  for  mathematical  modeling  of  the  dynamics 
for  the  optimal  controller  design  is  presented.  This  is  an  automated 
process  which  has  been  applied  to  multivariable  systems. 

System  performance  analysis  in  state  space  is  briefly  reviewed  in  Section  IV 
for  completeness.  The  steady-state  response  modeling  is  considered  in 
detail  since  it  is  required  for  the  ALDCS  design.  Some  of  the  analytical 
developments  reported  have  not  been  incorporated  into  the  existing  soft¬ 
ware  due  to  lack  of  resources. 

In  Section  V  active  control  synthesis  procedures  are  reviewed.  A 
description  of  performance  criteria  is  given  first.  This  is  followed  by 
the  control  configuration  for  the  required  steady-state  performance. 
Subsequently,  the  full  state  and  simplified  optimal  controller  design  steps 
are  described.  Finally,  a  demonstration  example  is  given  for  the  ALDCS 
design  using  the  C-5A  vehicle. 

In  Section  VI  the  effects  of  modeling  and  model  reduction  procedures  on 
control  system  performance  is  presented.  Two  vehicle  modeling  procedures. 


2 


i.e. ,  Air  Force  FLEXSTAB  and  Honeywell /GELAC,  are  reviewed  first 
for  the  C-5A  vehicle.  Subsequently,  the  theory  of  model  reduction 
procedures  is  briefly  described.  The  residualization  and  truncation 
methods  are  investigated  using  C-5A  open-loop  and  closed-loop  models. 

Conclusions  and  recommendations  are  given  in  Section  VII.  Both  future 
analytical  work  and  software  development  work  are  described. 

Computer  programs  which  implement  the  mathematical  analysis  and 
models  presented  in  this  volume  are  listed  in  Volume  II.  Complete 
documentation  of  KONPACT  is  beyond  the  scope  of  this  contract. 

In  Volume  III  users'  information  on  KONPACT  is  given.  The  input  cards 
are  fully  described  for  each  program.  A  brief  description  of  programs 
and  the  information  flow  in  KONPACT  are  also  presented  for  completeness. 
Demonstration  examples  are  included  to  guide  the  user  in  data  mechanics. 


SECTION  II 

ACTIVE  CONTROL  SYNTHESIS  APPROACH 

This  section  presents  the  overview  of  Honeywell's  work  on  the  optimal 
and  suboptimal  active  control  design  for  the  C-5A  aircraft  and  design 
software. 

First,  the  previous  work  of  Honeywell  on  the  subject  (Reference  4)  is 
briefly  presented  to  provide  background  information.  Subsequently  an 
overview  of  the  KONPACT  package-- Computer  Programs  for  Active 
Control  Technology,  developed  in  this  program  is  given. 

Finally,  the  results  of  the  Active  Lift  Distribution  Control  System  (ALDCS) 
using  the  FLEXSTAB/LSA  generated  vehicle  model  are  described. 

REVIEW  OF  DESIGN  PROCEDURE  FOR  THE  C-5A  ACTIVE 
LIFT  DISTRIBUTION  CONTROL  SYSTEM  (ALDCS) 


Honeywell  has  conducted  a  study  (Reference  4)  which  was  part  of  Lockheed 
Georgia’s  (GELAC)  program  to  increase  the  life  of  the  C-5A  airplane 
with  the  use  of  active  control.  The  design  goals  for  the  ALDCS  are  listed 
in  Table  1.  The  design  rules  were  to  meet  these  specifications  with  a 
single  accelerometer  in  each  wing  (see  Figure  1)  and  without  the  use  of 
notch  filters  to  remove  unwanted  flexure  modes  from  the  measurements. 
Four  flight  conditions  were  chosen  as  shown  in  Table  2. 


Wing  Panel  (2) 


Pitch  Rate  Gyro 
Accelerometer  — . 


Table  1.  ALDCS  Design  Goals 


No. 

Performance  Criterion 

Specification 

1 

Wing  root  bending  moment  due  to 
maneuvers  and  wind  gusts 

70%  of  free  A/C  values 

2 

Wing  root  torsion  due  to  wind  gusts 

Less  than  105%  free  A/C 

3 

Handling  qualities 

Same  as  A/C  with  existing 
SAS 

4 

Stick  force  per  g 

Same  as  A/C  with  existing 
SAS 

5 

Stability  margins 

10  db  gain,  60°  phase 

Table  2,  Flight  Condition  Definition 


Flight  Condition 

KONPACT 

CHECKCASE 

CRUISE 

412301 

CLIMB 

412502 

TRAFFIC 

412530 

CONTOUR 

412020 

F.C.  37 

Total  Weight  (lb) 

578,430 

698,400 

698,400 

529,500 

593, 154 

Mach  No. 

.74 

.448 

.228 

.333 

.533 

Altitude  (ft) 

30,000 

7,  500 

1,500 

0 

10,000 

Dynamic  Pressure 
(psf) 

240 

191 

73 

418 

292 

Airspeed  (fps)  Uq 

735 

468 

254 

594 

577 

Fuel  (lb) 

94,250 

214,500 

214, 500 

94.250 

159,750 

Cargo  (lb) 

158, 104 

160,000 

160, 000 

30.000 

110,000 

Center  of  Gravity 
(%  M.  A.  C. )  (c.g. ) 

30 

31 

30 

29.9 

31 
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The  resulting  system  was  called  Active  Lift  Distribution  Control  System 
(ALDCS).  It  was  developed  as  a  means  to  reduce  fatigue  damage  on  the 
wing  root  due  to  aircraft  flight  maneuvers  and  wind  gust  forces.  The 
control  system  developed  by  Honeywell  (Reference  4)  was  designed  to  reduce 
the  fatigue  damage  enough  to  double  the  actual  flight  life  of  the  aircraft. 

The  proposed  system  consisted  of  five  sensors  controlling  the  ailerons 
and  the  inboard  elevators.  An  accelerometer  was  placed  in  each  wing; 
an  accelerometer  and  a  pitch  rate  gyro  were  placed  in  the  fuselage,  and 
a  sensor  was  placed  on  the  pilots  pitch  control  column  (see  Figure  2). 

The  control  system  reduces  wing  bending  during  maneuvers  by  putting  an 
upward  bias  on  the  ailerons  proportional  to  the  g  load.  This  causes  the 
life  to  be  distributed  inboard  on  the  wing  and  requires  the  elevator  to  deflect 
more.  The  wing  bending  due  to  wind  gusts  and  maneuver  are  reduced  by 
sensing  these  forces  with  the  accelerometer  and  cancelling  them  with  the 
aerodynamic  surfaces.  The  aileron  most  effectively  damps  the  higher  fre¬ 
quency  bending  modes.  The  elevator  most  effectively  damps  the  low  fre¬ 
quency  bending  modes. 

Structural  mode  control  technology  has  been  developed  for  aircraft  in 
programs  with  the  XB-70,  B-52,  C-5A  and  YF-12A  (References  11  and 
18  through  21).  The  XB-70  program  demonstrated  the  effectiveness  of 
the  ILAF  (identical  location  of  accelerometer  and  force)  configuration 
where  the  sensor  and  control  force  producers  are  placed  near  one  another. 
The  control  surfaces  effectively  cancelled  out  sensed  aeroelastic  forces  due 
to  lower  local  acceleration. 
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Maneuver  Load  Control  and  Aileron 

Gust  Alleviation  Filter  ALDCS  Off  Actuator 


Figure  2.  ALDCS  Simplified  Controller  Configuration  for  Cruise  -  F.C.  412301 


The  B-52  and  C-5A  LAMS  programs  established  the  feasibility  of  using 
control  forces  to  reduce  fatigue  damage  ratio  in  air  turbulence.  The 
C-5A  LICS  program  showed  the  special  effectiveness  of  the  inboard  spoilers 
in  reducing  wing  root  fatigue.  The  YF-12A  study  extended  the  LAMS 
technology  by  examining  additional  control  force  producers.  A  canard 
vane  proved  to  be  the  most  effective  force  producer  for  reducing  structural 
mode  contributions  to  local  accelerations. 


In  the  C-5A  ALDCS  program  all  aspects  of  the  control  problem  were  considered. 
These  included  handling  qualities,  flutter  margins,  and  cost  effective  sensor 
configurations.  The  elevators  Droved  to  be  effective  force  producers  in 
combination  with  the  ailerons  for  reducing  wing  root  bending  moments.  In 
addition  to  low  frequency  bending  mode  suppression,  the  ailerons  were  up- 
rigged  for  steady-state  g  loads  to  distribute  the  maneuver's  lift  forces 
inboard  to  reduce  root  stress. 

Honeywell-GELAC  C-5A  Vehicle  Modeling 

Modeling  Procedure-One  of  the  most  difficult  and  crucial  tasks  in  the 
design  of  the  flexure  control  system  is  constructing  an  accurate  model  of 
the  aircraft  (Reference  5).  The  Lockheed  Georgia  Company  provided  the 
complex  model  of  the  C-5A  structural  dynamics  with  15  flexure  modes 
(Reference  4).  The  model  allowed  accurate  determination  of  flexure  mode 
states  over  the  aircraft  wing  as  a  function  of  aerodynamic  and  control  sur¬ 
face  forces.  The  provided  data  also  allowed  the  calculation  of  accelerometer 
outputs  for  any  sensor  location.  In  addition,  shear,  torsion,  and  bending 
moment  equations  were  provided  to  calculate  the  control  performance.  The 
above  data  were  supplied  for  the  four  critical  flight  conditions  presented  in 
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The  model  was  transformed  from  second  order  airframe  standard  form  to 
first  order  differential  equations  for  state  space  analysis.  Other  minor 
changes  were  made  to  the  model  to  simplify  the  state  space  analysis.  First 
order  actuator  models  were  added.  The  Von  Karman  wind  gust  model  was 
replaced  by  the  second  order  Dryden  model.  The  Wagner  dynamics  for  the 
flexure  modes  were  reduced  from  the  second  to  first  order  to  represent  un¬ 
steady  model.  The  phugoid  dynamics  were  removed  so  that  the  unstable  roots 
would  not  affect  the  design  calculations.  After  the  above  transformations, 
the  steady  and  unsteady  aircraft  models  contained  42  and  79  states, 
respectively. 

In  order  to  make  the  design  optimization  cost  effective,  the  design  model 
was  obtained  by  reducing  the  7th  through  15  flexure  modes  in  the  steady 
model  through  residualization.  By  this  process  the  highest  frequency  flexure 
modes  were  removed.  This  procedure  is  more  accurate  for  SS  maneuver 
load  calculations  than  modal  truncation. 

Specific  modeling  software  was  developed  to  perform  the  above  data  manipu¬ 
lations,  transformations,  and  reductions  (Reference  26).  We  note,  however, 
that,  because  of  the  way  the  unsteady  effects  (Wagner  dynamics)  were  modeled, 
the  steady  low  order  models  could  not  be  obtained  directly  from  the  high 
order  unsteady  model  either  by  truncation  or  by  the  residualization  processes. 

This  fact  was  overlooked  in  some  current  programs  at  Honeywell  using  the 
C-5A  high  order  model.  It  was  also  a  source  of  some  confusion  in  this  pro¬ 
gram  when  the  comparison  of  Honeywell/GELAC  and  FLEXSTAB/LSA  models 
were  made.  The  late  discovery  of  this  fact  prevented  reactivating  specific 
modeling  software  reported  in  Reference  26  to  obtain  accurate  reduced  order 
models  using  the  GELAC  data  for  this  study. 


Coordinate  Systems--The  rigid  body  dynamics  of  the  Honeywell  study  are 
described  in  the  body  fixed  axis  systems.  The  origin  of  this  system  is  at 
the  aircraft's  center  of  gravity  (c.  g. )  and  moves  with  the  c.  g.  along  the 
flight  path.  (See  Reference  4  for  additional  details.) 

Sending  moments  and  torsional  moments  provided  by  Lockheed  are 
described  at  the  five  wing  station  local  coordinate  systems  shown  in 
Figure  3.  Polarity  is  established  by  the  "  right-hand"  sign  convention. 

Bending  moments  are  about  an  axis  perpendicular  to  the  elastic  axis 
(+  bending  moment  produces  right  wina-tip  up) .  Torsional  moments  are 
about  an  axis  parallel  to  the  elastic  axis  (  +  torsion  produces  leading- 
edge  up). 

Active  Control  Using  Quadratic  Optimization 

The  quadratic  optimization  design  methods  were  carried  out  on  a  six  mode 
model  without  Wagner  dynamics.  The  design  model  required  28  states  with 
various  compensator  states  included. 

Quadratic  optimization  is  a  numerical  technique.  The  technique  is  especially 
suited  to  complex  problems  with  many  interactive  loops  and  many  conflicting 
performance  criteria.  In  such  problems,  intuition  is  confused.  The 
numerical  technique  provides  a  systematic  procedure  for  making  complex 
tradeoffs.  By  simply  adding  all  the  performance  parameters  to  a  response 
vector  and  weighing  their  relative  importance,  the  method  gets  near  the 
desired  result  very  quickly. 

For  quadratic  optimization  the  model  and  its  responses  are  reduced  to 
algebraic  quantities.  The  rms  responses  are  described  by  the  coefficients 
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Bi  =  Bending  Moment 
Ti  =  Torsional  Moment 


Figure  3.  Equivalent  Honeywell /GELAC  C-5A  Load  Axis  Systems 


Winq  Station  (In) 


Bi  =  Bending  Moment 
Ti  -  Torsional  Moment 


Figure  4.  Equivalent  FLEXSTAB  C-5A  Load  Axis  Systems 


of  the  differential  equations  of  the  system.  These  responses  can  be  summed 
in  a  cost  function  with  variable  weights.  The  feedback  gains  are  chosen  to 
minimize  the  sum  of  the  responses.  The  design  problem  is  to  pick  the 
proper  responses  for  the  cost  function  with  the  proper  weight  ratio  for  the 
performance  tradeoff.  Initially,  feedback  gains  are  chosen  for  all  the 
states  (dynamic  elements)  of  the  system.  Since  only  a  limited  number  of 
states  can  be  fed  back  through  the  sensors,  the  gains  have  to  be  adjusted 
later  for  that  condition  in  order  to  make  the  system  practical. 


Cost  Function 


The  cost  function  for  the  quadratic  design  was  made  up  of  flexure  mode 
rates,  bending  and  torsion  moments  at  several  wing  stations,  control 
surface  position  and  velocity,  and  pitch  response  model  error.  All  of 
the  above  responses  were  minimized  in  accordance  with  selected  weighting 
ratios.  The  pitch  response  of  the  C-5A  with  the  stability  augmentation 
system  (SAS)  was  used  as  a  model  for  handling  qualities. 

Design  Procedure 


The  quadratic  optimization  of  the  above  cost  function  produced  a  full  state 
feedback  design.  The  design  was  then  "practicalized"  by  adjusting  the 
feedback  gains  to  be  compatible  with  the  actual  sensors.  This  process  is 
accomplished  with  a  numerical  search  algorithm  which  seeks  to  adjust  the 
feedback  gains  with  the  least  amount  of  increase  in  the  cost  function. 
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In  this  design  procedure,  for  each  response  weights,  the  DIAK  program 
generates  full  state  feedback.  Successive  response  weights  are  chosen 
until  performance  specifications  are  exceeded.  Then  the  design  is  practicalized 
using  the  FFOC  program.  If  the  practicalized  design  meets  performance 
specifications,  performance  plots  are  calculated  using  the  analysis  model. 

If  not,  a  new  full  state  design  is  calculated.  If  the  frequency  response  plot 
shows  insufficient  stability  margin,  slight  gain  or  filter  changes  are  made. 

Then  the  rrns  responses  and  transient  responses  are  rechecked. 

Performance  Evaluation 


During  the  design  process,  the  performance  of  the  system  was  evaluated 
by  examining  the  rms  values  of  the  wing  root  bending  and  torsion  moments. 
The  handling  qualities  were  evaluated  by  comparing  the  ALDCS  design  tran¬ 
sient  response  to  a  step  pilot  command  with  similar  SAS  aircraft  responses. 
In  addition,  the  damping  ratios  and  the  natural  frequency  of  the  short  period 
roots  are  compared  with  that  of  the  SAS  aircraft.  The  stability  of  the  ALDCS 
system  is  qualitatively  evaluated  by  looking  at  the  damping  ratios  of  all  of 
the  roots  of  the  characteristic  equation.  A  quantitative  measure  of  stability 
margin  is  obtained  from  loop  breaking  frequency  response  plots.  One  sensor 
or  actuator  loop  is  opened  at  a  time  to  determine  the  gain  margin  for  each 
loop. 

ACTIVE  CONTROL  TECHNOLOGY  COMPUTER  PROGRAM 

4  computer  program  (KONPACT)  was  generated  to  facilitate  dynamic 
modeling,  optimal  and  suboptimal  controller  synthesis,  and  performance 
analysis  of  vehicles  with  lightly  damped  modes. 
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This  program  utilizes  advanced  computational  techniques  to  perform  system 
modeling,  modern  control  synthesis,  analysis,  and  design  of  automatic  con¬ 
trol  systems.  Figure  5  shows  its  functional  block  diagram.  It  interfaces 
with  the  FLEXSTAB/LSA  program  for  vehicle  description  as  well  as  per¬ 
formance  evaluation  of  the  optimal  closed  loop  system.  Figure  6  shows  the 
LSA/KONPACT  interface  and  also  shows  a  proposed  version  of  an  advanced 
FLEXSTAB  program.  KCNPACT  is  also  used  as  a  stand  alone  program  with 
externally  input  vehicle  descriptions.  Specifically,  KONPACT  integrates 
the  relevant  computer  programs  in  aerodynamics,  and  structures  (FLEXSTAB/ 
LSA  level  2.  01. 00)  with  programs  in  modern  control  theory  (DIAK,  FFOC) 
into  a  single  interdisciplinary  design  tool.  This  integration  is  shown  in  detail 
in  Figure  7.  Its  variable  dimensioning  feature  reduces  the  workload  on  the 
interface  data  mechanics  for  design  and  analysis  of  large  scale  systems. 

In  the  following  a  brief  description  of  KONFhCT  is  given. 

Overlay  Organization 

KONPACT  consists  of  two  programs,  namely,  a  modeling  program 
(KONPACT-1)  and  a  design  program  (KONPACT-2).  KONPACT-1  interfaces 
with  FLEXSTAB  through  the  LSA  program  to  obtain  the  vehicle  model  and 
augments  the  specified  dynamics  to  obtain  the  state  space  description 
(quadruple  data)  of  the  flight  control,  system.  These  data  are  utilized  by 
KONPACT-2  which  contains  the  subprograms  DIAK  and  FFOC  (documented 
in  Reference  1)  in  the  design  of  the  optimal  feedback  gains.  Also  KONPACT-2 
interfaces  with  FLEXSTAB  through  the  LSA  program  to  evaluate  performances 
of  the  above  designed  optimal  flight  control  system. 
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Figure  7.  KONPACT  Program  Organization 
(Dotted  lines  show  future  work) 


Table  3  provides  a  brief  description  of  programs  KONPACT-1  and  KONPACT-2 
and  their  subprograms.  The  detailed  interface  between  KONPACT  and  the 
LSA  program  is  illustrated  in  Figure  8.  The  overlay  structure  of  KONPACT 
is  discussed  in  References  51. 

Information  Flow 

The  normal  sequence  for  obtaining  an  overall  state  space  model  of  a  flight 
control  system  using  the  modeling  program  (KONPACT-1)  is  as  follows: 

•  The  vehicle  model  is  obtained  by  using  either  subprogram  STAMK1 
for  LSA  data  or  subprogram  STAMK4  for  other  types  of  vehicle 
data. 

•  The  actuator,  sensor,  controller,  implicit  and  explicit  models 
are  obtained  by  using  either  subprogram  STAMK2  with  transfer 
function  input  data  or  subprogram  STAMK3  with  quadruple  input 
data. 

•  The  subsystems  defined  above  are  combined  to  get  an  overall 
system  by  using  subprogram  STAMK3  with  interconnection  input 
data. 

•  The  overall  system  model  is  conditioned  (modified)  by  scaling 
and/or  shuffling  and/or  truncating  and/or  residualizing  the 
variables  using  the  CONDK  program.  This  program  also  develops 
the  rate  of  change  of  response  variables  v/hen  required. 

The  normal  sequence  for  designing  optimal  feedback  controllers  and  evaluating 
the  performance  of  the  resulting  system  using  the  design  program  KONPACT-2 
is  as  follows: 


Table  3.  KONPACT  Program  Descriptions 


PROGRAM 

SUBPROGRAM 

DESCRIPTION 

KONPACT- 1 

State  space  modeling  program 

STAMK.1 

Obtains  state  space  model  from  LSA 
simulator  deck  data 

STAMK2 

Obtains  state  space  model  from  transfer 
function  data 

STAMK3 

Obtains  state  space  model  from 
quadruple  data  and  interconnection  data 

STAMK4 

Obtains  state  space  model  from 
simulation  equations  (User  Written) 

CONDK 

Modifies  the  state  space  model  by 
scaling,  shuffling,  truncating  and 
residualizing  the  system  variables 

KONPACT-2 

Optimal  design  program 

DATAK 

Prepares  data  for  DIAK,  FFOC  and 

LSA  programs 

DIAK 

Designs  full  state  feedback  optimal 
controllers 

FFOC 

Designs  reduced  state  (practical.) 
feedback  optimal  controllers 
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Figure  8.  Interface  Between  LSA  and  KONPACT  Programs 


•  Punching  Simulator 
Interface  Deck  for 
Vehicle 


•  Obtaining  Frequency 
Domain  Representation 
of  Vehicle 


Pole  and  Zero  Location 
Evaluation 

Frequency  Response 
Evaluation 

Power  and  Power  Spectral 
Density  Evaluation 


•  Full  state  feedback  control  gains  are  obtained  by  using  the  DLA.K 
program  and  by  varying  the  quadratic  weights  until  performance 
design  requirements  are  satisfied. 

•  The  resulting  full  state  feedback  control  gains  are  reduced  to  gains 
only  on  specified  measurements  by  using  the  FFOC  subprograms. 

•  The  performance  of  the  resulting  closed  lcop  system  is  evaluated 
using  the  LSA  program. 

•  The  above  steps  are  repeated  until  a  satisfactory  design  is  obtained. 

Table  4  describes  all  the  data  tapes  used  in  the  KONPACT-1  and  KONPACT-2 
programs.  The  state  space  model  data  (quadruple  data)  and  the  name  list 


Table  4.  KONPACT  Data  Tapes 


'  TAPE 
NAME 

GENERATING 

PROGRAM 

BENEFITING 

PROGRAM(S) 

VDATA 

Simulator  Interface  data  in  the 
form  of  card  images 

LSA 

KONPACT-1 

QDATA 

Quadruple  (A,B,C,D)  or  state 
variable  representation  data 

KONPACT-1 

KONPACT-1 

KONPACT-2 

NDATA 

Name  list  data  of  the  state 
variable  representation 

KONPACT-1 

KONPACT-1 

DDATA 

Full  state  feedback  gain  data 
in  the  form  of  card  images 

KONPACT-2 

KONPACT-2 

FDATA 

Reduced  feedback  gain  data 
in  the  form  of  card  images 

KONPACT-2 

KONPACT-2 

SDSTP 

Frequency  domain  representa¬ 
tion  of  quadruple  data 

KONPACT-2 

LSA 

data  are  written  on  tapes  QDATA  and  NDATA,  respectively.  The  vehicle 
data  (simulator  deck  data)  are  written  on  tape  VDATA.  The  feedback  gain 
data  from  DTAK  and  FFOC  are  written  on  tapes  DDATA  and  FDATA, 
respectively.  The  overall  system  data  in  frequency  representation  form 
are  written  on  tape  SDSTP  for  use  by  the  LSA  program. 

The  DATAK  subprogram  is  used  in  preparing  data  tapes  for  DIAK,  FFOC, 
and  LSA. 

Variable  Dimensioning 

Dynamic  data  storage  variable  dimensioning  techniques  (Reference  7) 
are  used  for  efficient  data  storage.  This  technique  also  facilitates  changing 
the  amount  of  allocated  (required)  storage  space  by  a  data  card  input. 

In  KONPACT  the  subprogram  arrays,  whose  size  depends  on  the  maximum 
system  dimension  inputs,  are  stored  in  scratch  storage  blocks  using 
variable  entry  points.  In  the  subprograms  the  arrays  are  dimensioned 
with  integer  variables.  These  "variable  DIMENSION  statements"  remain 
unchanged  although  the  amount  of  required  data  storage  is  altered.  The 
maximum  si  :e  of  the  scratch  storage  blocks  is  specified,  in  a  "fixed 
DIMENSION  statement,  "  in  the  main  program. 

The  size  of  storage  actually  needed  by  the  arrays  varies  depending  on  the 
maximum  system  dimension  inputs.  Thus,  if  the  maximum  size  a  user 
allows  for  his  program  changes,  there  is  only  the  "fixed  DIMENSION 
statements"  in  the  main  program  to  be  changed.  Changing  the  main  program 
of  KONPACT- 1  is  done  by  a  precompiler.  The  user  provides  the  new 
maximum  system  dimensions  by  data  cards.  Updating  and  running  with  the 


updated  main  program  are  done  with  control  cards  in  a  single  run.  For 
more  details  on  variable  dimensioning  the  reader  is  referred  to  Volume  II 
(Reference  3). 

ALDCS  DESIGN  WITH  FLEXSTAB  C-5A  VEHICLE  DATA 

Two  of  the  contract  objectives  were  1}  to  check  out  KONPACT  by  repeating 
the  past  Honeywell  design  work  (Reference  4)  using  the  FLEXSTAB  generated 
C--5A  vehicle  data  and  2)  to  compare  and  correlate  the  resulting  FLEXSTAB 
controller  to  a  Reference  4  ALDCS. 

Modeling  Procedures 

The  Air  Force  supplied  C-5A  data  for  the  cruise  flight  condition  in  the  form 
of  cards  (simulator  data  deck).  The  FLEXSTAB  simulator  data  were  con¬ 
verted  to  state  format  and  augmented  by  the  wind  and  wind  distribution  states, 
actuators,  and  controller  dynamics  as  shown  in  Figure  9.  Subsequently, 
the  augmented  data  were  scaled  and  shuffled  to  correspond  with  Honeywell/ 
GELAC  data. 

Controller  Design 

After  having  obtained  agreement  between  the  two  sets  of  data,  a  full  state 
design  was  obtained  using  the  weights  of  Reference  4  and  the  DIAK  program 
in  KONPACT  software. 

Subsequently,  the  gains  on  the  measured  variables  were  retained  and  other 
gains  were  reduced  to  zero  using  the  FFOC  program  in  KONPACT  software. 


Figure  9.  General  Overall  System  Block  Diagram 


Finally  the  pilot  gain  the  reduced  control  law  was  adjusted  to  maintain  the 
same  stick  gradient  per  g  as  that  of  SAS.  The  closed-loop  model  was  pro¬ 
duced  using  KONPACT  software  with  the  final  ALDCS  controller  to  evaluate 
several  performances. 

Time  history  plots  were  made  using  DIAK.  Power  spectral  density  plots 
were  made  using  the  FLEXSTAB/LSA  program.  These  results  are  presented 
in  detail  in  Section  V. 

This  repeat  design  process  thoroughly  checked  out  the  total  KONPACT  system 
with  respect  to  communication  with  its  internal  programs  as  well  as  with  the 
FLEXSTAB/LSA  program  system. 

Coordinate  Systems 

The  FLEXSTAB  rigid  body  dynamics  are  described  in  the  body  fixed  axis 
system.  The  origin  of  this  system  is  at  the  aircraft's  c.  g.  and  moves  with 
the  c.g.  along  the  flight  path.  (See  Reference  2  for  additional  details.) 

Bending  moments  end  torsional  moments  are  described  at  the  five  wing 
station  local  coordinate  systems  shown  in  Figure  4.  The  polarity,  orientation, 
and  location  of  these  load  axis  systems  are  equivalent  to  the  Honeywell/GELAC 
data  with  one  exception:  the  FLEXSTAB  inboard  (w.s.  =  120  in)  system  is 
rotated  18.6°  right  wing  tip  forward.  Due  to  the  late  discovery  of  this 
fact,  all  design  performed  with  the  FLEXSTAB  model  includes  this  discrepancy. 
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SECTION  III 


DYNAMIC  MODELING  FOR  OPTIMAL 
CONTROLLER  DESIGN 


In  this  section  we  briefly  present  procedures  on  the  automatic  modeling 
of  interconnected  dynamic  systems  for  optimal  control  synthesis  and 
analysis.  State  variable  form  is  selected  to  characterize  the  dynamic 
elements  of  a  system.  This  is  necessary  for  optimal  control  synthesis. 

Also  a  uniformity  in  the  model  form  (irrespective  of  the  size  or  the  internal 
structure  of  subsystems)  facilitates  the  analysis  during  the  evaluation  of 
various  performance  measures. 

In  the  state  space  representation  a  dynamic  model  is  characterized  by 
four  matrices  (A,B,C,  D)  for  the  continuous  model  and  (F,  G,  H,E)  for  the 
discrete  model.  These  four  matrices  are  referred  to  as  a  system  quadruple. 
First  we  present  an  approach  to  develop  a  system  quadruple  from  the 
physical  equations  describing  the  dynamics  of  a  system.  Converting  the 
FLEXSTAB/LSA  simulator  data  matrices  to  quadruple  format  is  treated 
here.  This  is  followed  by  an  approach  for  modeling  with  transfer  function 
inputs.  This  facilitates  the  augmentation  of  vehicle  dynamics  with  the 
actuator,  filtered  measurements,  and  controller  dynamics.  Subsequently 
modeling  of  transport  delays  in  the  vehicle  aerodynamics  is  treated  using 
the  Pade*  approximation  and  developing  the  corresponding  quadruple. 

This  is  followed  by  the  overall  system  modeling  using  the  dynamic  blocks 
and  the  interconnections  between  them.  The  response  rate  modeling  is 
treated  next.  Successful  optimal  control  synthesis  requires  careful 


construction  of  design  responses  to  force  the  desired  performance.  The 
desired  performance  in  general  contain  responses  and  response  rates. 

Finally  we  present  briefly  the  condition1  ng  of  a  model  for  optimal  control 
design.  Model  conditioning  includes  reduction,  shuffling,  and  scaling 
operations.  Model  reduction  can  be  achieved  in  several  ways.  Here  we 
present  only  the  residualization  and  truncation  procedures. 

Reordering  of  states  is  called  shuffling.  This  operation  is  needed  to  bring 
two  different  model  variables  into  a  common  base  for  comparison  as  well 
as  to  reduce  the  system  order.  The  scaling  operation  also  facilitates 
the  data  comparison  and  overall  system  modeling  by  bringing  subsystem 
units  into  a  common  base. 

DEVELOPMENT  OF  THE  LINEAR  SYSTEM  MATRICES 
FROM  THE  SIMULATION  EQUATIONS 

In  general,  the  simulation  equations  of  the  system  take  the  following  form: 


X  =  f(x,  y. 

X,  u) 

(1) 

y  =  g(x,  y. 

X,  u) 

(2) 

it 

X* 

«• 

x,  u) 

(3) 

where 

x  =  n  x  1  vector  of  the  output  of  integrators 
y  =  n  x  1  vector  of  the  output  of  summing  points 


r  =  x  1  vector  of  the  system  variables  of  interest  (response  outputs) 
u  =  n^  x  1  vector  of  the  external  inputs 


The  functions  f,  g  and  h  are  usually  nonlinear.  For  the  linear  analysis 
they  can  be  linearized  about  a  given  operating  point.  In  the  following, 
we  shall  assume  that  the  simulation  equations  represent  the  linearized 
model.  In  this  case.  Equations  (1),  (2)  and  (3)  can  be  put  in  the  following 
form: 

x  =  F.  x  +  F  y  +  F  x  +  F  u 
x  y  J  x  u 

y  =  G.  x+G  y  +  G  x+G  u 
x  y  x  u 

r  =  H.  x+H  y  +  H  x+H  u 
x  y  J  x  u 

and  this  set  of  equations  can  be  reduced  to  the  following  standard  form 
by  algebraic  operations 

x  =  Ax  +  Bu 
r  =  Cx  +  Du 

On  the  surface,  this  task  appears  to  be  very  simple  to  carry  out  with 
paper  and  pencil.  However,  for  large  systems  the  writing  of  simulation 
equations  in  the  format  given  in  Equations  (4),  (5),  and  (6)  is  prone  to 
human  error  and  should  be  avoided. 


(7) 

(8) 


(4) 

(5) 

(6) 


In  the  following,  we  present  an  algorithm  which  automates  the  transition 
from  the  physical  equations  (analog  simulation  equations)  to  the  state 
variable  representation  given  by  Equations  (7)  and  (8). 
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Let  us  define  two  vectors  as  follows. 


v  =  col  (x,  y,  r)  (9) 

w  =  col  (x,  y,  x,  u)  (10) 

Obviously,  Equations  (4),  (5)  and  (6)  can  be  written  as 

v  =  F(w)  (11) 


The  matrix  coefficients  given  in  Equations  (4),  (5)  and  (6)  are  then  obtained 
by  first  finding 

bF 
bw 

x  p 

and  then  properly  partitioning  it.  This  term  —  is  called  the  simulation 

bw 

matrix.  The  sizes  of  its  rows  and  columns  are  given  respectively  by 


n  =  n  +  n  +  n 
x  y  r 


(12) 


m  =  2n  +  n  +  n 
x  y  u 


(13) 


The  coefficient  matrices  obtained  by  partitioning  the  simulation  matrix 
is  indicated  in  Figure  10. 
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The  column  vectors  i  =  1,  2,  . ..  ,  m  are  obtained  simply  by  setting 


w.  =  1 

l 


(14) 


y  0,  j  =  1,  2 .  m,  j  H 

and  evaluating  (11).  This  yields  the  coefficient  matrices. 


In  the  sequel,  the  algebraic  reduction  process  will  be  described.  First, 
Equations  (4)  and  (5)  are  written -It?  the  following  form: 


r a  -  f  j ! 

-F 

lk\ 

'  F 

F 

X 

y 

= 

X 

u 

-G. 

(I  -  G  ) 

1  y 

G 

G 

X 

y 

\  i 

X 

u  _ 

(15) 


Then 


/.\ 

x 


y 

\  / 


is  obtained  in  terms  of  x  and  u  by  solving  Equation  (15). 
Then  r  is  obtained  in  terms  of  x  and  u  by  substituting  (15)  into  (16): 


r  =  (H  |  H  ) 
x  y 


x 


LyJ 


+  (H  !  H  ) 
x  1  u 


(16) 


The  subroutine  which  implements  this  algorithm  is  called  STAMK. 
Implementation  of  the  Simulation  Equations 

The  physical  (simulation)  equations  describing  the  system  dynamics 
(Equations  (4),  (5)  and  (6))  are  implemented  in  subroutines  SIMK1  and 
SIMK2. 
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The  physical  equations  describing  the  vehicle  dynamics  in  the  FLEXSTAB 
system  (Reference  2)  are  shown  in  Volume  III,  Figure  9  (Reference  51). 
These  equations  are  implemented  in  subroutine  SIMK1.  The  coefficient 
matrices  are  read  from  the  simulator  interface  deck  in  the  beginning  of  the 
above  subroutine. 

MODELING  WITH  TRANSFER  FUNCTION  INPUT 

In  the  following  we  present  an  approach  to  carry  out  system  modeling  by 
software  with  transfer  function  input.  The  approach  consists  of  two  parts: 
1)  obtaining  the  corresponding  quadruple  for  each  transfer  function  block, 
and  2)  combining  the  blocks  using  the  connection  equations  and  obtaining 
the  system  quadruple.  In  the  following  we  discuss  each  in  that  order. 

Consider  a  system  characterized  by  its  output/input  relation: 


R(s)  _  w/  . 
U(s)  H(s) 


b  sn  +  b  1sn_1  +  ...  +  h  s  +b 
n _ n-1  1  o 


n  ,  n-1 

ans  +  Is  +  ...  +  a.s  +  a 
n  n-1  1  o 


>  V0 


(17) 


There  are  many  ways  of  realizing  this  transfer  function.  (See  Reference 
31  for  major  realization  forms. )  In  the  following  we  shall  develop  input 
Frobenius  form  realization  and  obtain  corresponding  quadruple  in  para¬ 
metric  form  for  software  implementation. 


The  long  division  of  Equation  (17)  yields 


The  transfer  function  coefficients  in  Equation  (17)  form  a  2  x  n  array  as 


indicated  below. 


H(s)  = 


.  L1  _  J 
a  1 
n  i 


b  b 

a,  i  a 

-  1  i  o. 


Equations  (20)  and  (21)  form  an  algorithm  for  obtaining  the  quadruple  of 
th 

an  n  order  transfer  function.  Subroutine  TRANSK  implements  this 
algorithm. 


To  develop  the  system  quadruple,  one  must  combine  the  block  quadruples 
obtained  as  described  above  using  the  connection  relations.  To  demonstrate 
the  approach  taken,  consider  a  block  diagram  of  a  system  containing  three 
transfer  function  blocks  as  shown  in  Figure  12. 


Ml)  |<X’)  Ir.O) 


uf(3)  I  (<3> 


Figure  12.  Block  Diagram  of  a  System  Containing  Three 
Transfer  Function  Blocks 
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Each  block  is  identified  by  four  quantities:  1)  a  block  number,  2)  HS  array 
representing  the  transfer  function  data,  3)  state  number,  and  4)  output- 
input  pair.  We  note  that  the  inputs  and  outputs  (i.e.,  u(l),  u(2)  and  r(l)) 
external  to  the  box  are  unsubscripted  variables,  whereas  inside  the  box 
they  are  subscripted  with  i  denoting  that  they  are  internal  variables. 


With  these  definitions  the  simulation  equations  corresponding  to  the  system 
shown  in  Figure  12  can  be  written  as  follows. 


Xf  =  AjXj  +  B1u.(l) 


x^  =  A2x2  +  B2u.<2)  )  Dynamics 


X3  =  A3X3  +  B3Ui^3) 


where  HSi  is  defined  by  (A.  B.  C.  D.). 

riU)  =  C1X1  +  Diui(1) 

r.(2)  =  C_x0  +  D  u.(2)>  Internal  outputs 
i  2  2  2  1  [ 


ri(3)  =  C3x3  +  D3u.(3^ 


(22) 


(23) 


Ui  =  p  ri  + 


r  =  R  r.  +  Su 


Internal  inputs 


External  output 


(24) 

(25) 


J 
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The  set  of  equations  given  above  are  implemented  in  subroutine  SIMKT. 
The  quadruple  (P,  Q,  R,  S)  appearing  in  Equations  (24)  and  (25)  are  called 
the  connection  quadruple.  For  the  system  shown  in  Figure  12  their  values 
are  given  below. 


0 

1 

— * 

1 

0 

0 

p  = 

-1 

0 

0 

9 

Q  = 

1 

0 

0 

0 

°. 

0 

1 

R  = 

(1 

0 

0) 

» 

s  = 

(0 

0) 

The  system  quadruple  is  obtained  via  STAMK  as  described  previously. 


MODELING  FOR  TRANSPORT  DELAYS 


The  transport  delays  resulting  from  the  gust  penetration  effect  in  the 
development  of  aerodynamic  forces  and  moments  are  represented  by 

H(s)  =  eSTi  i  =  1,2,3  (27) 

where  T^  is  the  time  delay  at  the  i  gust  input  station.  This  irrational 
transfer  function  in  s -plane  may  be  converted  to  a  rational  transfer  function 


H.(s) 


N(s) 

D(s) 


(28) 


in  different  ways.  Here  we  adopt  the  Pade'  approximation  (Reference  30) 
for  this  conversion,  and  Table  5  provides  the  numerator  and  denominator 
polynonial  coefficients  for  a  specified  degree  of  these  polynomials.  The 
Pade^  table  is  inplemented  in  subroutine  DFN. 
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OVERALL  SYSTEM  QUADRUPLE 


To  develop  the  overall  system  quadruple  ,  one  must  combine  the  subsystem 
quadruple  using  the  interconnection  relations.  The  approach  is  similar 
to  the  previous  case.  The  main  difference  in  this  case  is  that  quadruples 
for  the  subsystem  are  part  of  the  input  data,  and  each  subsystem  may 
have  multiple  inputs  and  outputs.  Consider  the  block  diagram  of  the  system 
consisting  of  two  subsystems  as  shown  in  Figure  13,  Each  subsystem  is 
identified  by  four  quantities:  1)  system  number,  2)  quadruple  data  for  the 
system,  3)  system  states,  4)  system  outputs  and  inputs.  Again  the  inputs 
and  outputs  of  the  overall  system  are  unsubscripted  variables,  whereas 
inside  the  box  they  are  subscripted  with  i  denoting  that  they  are  internal 
variables. 


Figure  13.  Overall  System  Block  Diagram 
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The  simulation  equations  corresponding  to  the  system  can  be  written  as 
follows : 


x„  =  A„x,  +  B.u.. 
1  11  1  ll 


*2  =  A2X2  +  B2Ui2 


-s 


)  Dynamics 


J 


r. . 

=  C,x, 

+  D„u., 

ll 

1  1 

1  il 

r._ 

=  C„x„ 

+  D„u.„ 

i2 

2  2 

2  i2 

>  Internal  Outputs 


(29) 


(30) 


where  System  i  is  defined  by  quadruple  (A.B.C^D.). 


Uil  “  P11  ril  +  P12  ri2  +  V 


"N 


Ui2  =  P12  ril  +  P22  ri2  +  Q2U 


>  Internal  Inputs 


r  =  R1  ril  +  R2  ri2  +  Su 


} 


J 


External  Output 


(31) 


(32) 


The  above  set  of  equations  are  implemented  in  subroutine  SIMK.  The  set 

of  matrices  {P. Q.,  R.,  S}  are  called  the  interconnection  quadruples. 
x3  ^  3 

The  combined  system  quadruple  is  obtained  via  STAMK  as  before. 


MODELING  FOR  RATE  RESPONSES 


The  original  system  is  described  by 

x  =  Ax  +  Bu 
r  =  Cx  +  Du 

From  the  response  set  specification  the  program  computes 

r  =  C  x  +  D  u 
s  s  s 


(33) 

(34) 


(35) 


Then,  Equations  (33)  and  (35)  form  the  description  of  the  new  system  with 

specified  responses.  The  elements  of  C  and  D  matrices  are  constructed 

s  s 

as  shown  in  Table  6. 

As  can  be  seen  from  this  table,  unless  D,  is  a  null  matrix,  the  input  space 

J 

should  be  extended  to  include  u.  When  the  implicit  model  error  rate  response 
is  specified  in  the  response  set,  the  program  computes  it  using  the  following 
algorithm. 

Implicit  Model  Following  Error  Response 

In  general  the  description  of  the  original  system  (33)  and  (34)  contains  the 
implicit  model  as  follows: 


X 

m 


41 


A 

0 

X 

B 

P 

P 

+ 

P 

0 

A 

X 

B 

_ 

m_ 

m_ 

L  m 

u 


(36) 


Construction  of  Equation 


r 

P 

"cp 

X 

p 

1 - 

a 

+ 

r 

e 

C1  C2_ 

X 

m 

—  _ 

0 

(37) 


where  xp  is  the  plant  state  vector  and  *m  the  implicit  model  state  vector  and 

e  is  the  error  vector  in  the  responses  between  plant  and  implicit  model.  If 

the  specified  response  set  includes  implicit  model  error  rates,  we  obtain 
from  (36)  and  (37) 


'*  =  ClVp  +  C2Amxm  +  (ClBp  +  C2Bm>  « 


It  is  assumed  that  the  implicit  model  following 
approximately 


(38) 

error  is  small.  This  yields 


=  ClXp  +  C2Xm  =  0 


(39) 


or 


-1 


x  =  -C  C  x 

m  2  Vp 


(40) 


Substituting  this  in  (38)  we  get 

re  =  <ClAp  "  C2AmC2  1<:i)xp  +  <C2Bp  +  C^)  u 


(41) 


Finally  the  overall  system  description  with  implicit  model  following 
rates  is  given  by 


error 


x  =  A  x  +  B  u 
P  P  P  P 


(42) 
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Reduction  of  the  size  of  a  dynamic  equation  of  a  system  can  be  achieved  in 
several  ways  depending  upon  the  reduction  criteria. 


In  the  following  we  present  two  approaches  for  reduction.  Table  7  shows 
these  and  criteria  for  reduction. 


Mam 


r 

C, 

“  “ 

D 

P 

1 

x  + 

P 

P 

• 

r 

c 

D 

e 

e 

e 

_  _j 

—  — 

—  — 

(43) 


where 


D  =  (C.B  +  C0B  ) 
c  1  p  2  m 


(44) 

(45) 


If  the  specified  response  set  does  not  include  implicit  model  following 
error  rates,  then  the  implicit  model  states  are  directly  truncated  to  obtain 


(46) 

(47) 


MODEL  CONDITIONING  (REDUCTION,  SHUFFLING  AND  SCALING) 
Reduction  of  Order 


Table  7.  Reduction  Criteria  ( x 2  =  collection  of  states  to  be  reduced) 


Reduction  Procedure 

Criteria  for  Reduction 

Truncation 

x2  =  °,  x2  =  ° 

Residualization 

X2  =  x2ss’  *2  =  0 

Now  consider  the  system  described  by 
x  =  Ax  +  Bu 
r  =  Cx  +  Du 


(48) 


where  x  and  r  are  the  state  and  response  vectors  of  appropriate  sizes. 
Partition  the  state  vector  x  =  col  (x^,x2)  where  x2  represents  the  collection 
of  state  components  to  be  reduced.  Equation  (53)  can  be  written  as 


• 

A„  „ 

A„  _ 

x. 

B. 

1 

= 

11 

12 

1 

+ 

1 

_  X2  _ 

_A21 

A22. 

_X2. 

-B2_ 

X1 


(49) 


r  =  (C,  C0) 


+  Du 


The  truncated  model  is  obtained  by  setting  =  0  and  x 2  =  0.  The  reduced 


system  is  described  by  (~  denotes  reduced  system  variables) 


*i  ■  An  xi +  Biu 


(50) 


r  =  Xj  +  D  u 


Residualized  model  is  obtained  by  setting  =  0  and  x^  =  x2ss* 
This  implies 


0  ~  A21  X1  +  A22  x2  +  B2  U 


(51) 


Assuming  that  is  a  stable  matrix  we  get 


X2  =  “A22  ^A21  X1  +B2u) 


(52) 


The  substitution  of  (52)into  (49)'  yields  the  residualized  system  model  as 
follows : 


“  (Aji  "  Ai2A22  A21^  X1  +  ^B1  ~  A12  A22  B2^U 


(53) 


?  -  <Cj  -  C2A-1  A21)  +  (D  -  C2A22  B2)  u 


-1 


x2  =  A21  X1  +  A22  X2  +  B2U 


The  smaller  the  time  constants  associated  with  states  x2  (compared  with 


the  time  constants  associated  with  states  x^  the  closer  the  agreement  is 
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ffHHi 


between  the  time  responses  of  the  residualized  system  (53)  and  the  original 
system  (49)  to  step  inputs. 

Define 

J-V/ 

x  =  x  +  e  (54) 

/v 

where  x,  x  and  e  denote  actual  state  response,  approximate  state  response 
and  error  response,  respectively. 

Using  Equations  (49)  and  (53)  in  (54)  one  obtains  the  differential  equation 
of  the  residualization  error  as  follows 

(55) 

(56) 

with  the  initial  conditions  e^(0)  ~  0,  ££(0)  £  0  . 

Equation  (55)  shows  that  the  residualization  error  is  driven  by  the  input 
u  as  well  as  actual  system  states  x  and  x  . 

X  z 


e„  = 


Anei 


+  (A 


11 


-  An)  xx  +  A12x2  +  (B1  -  B±)  u 


e„  = 


A2 1 e  1 


+  A22e2 


The  integral  of  the  square  error  can  be  computed  by  augmenting  (55)  to  (40) 
and  using  the  Lyapunov  equation  (Reference  31).  The  resulting  performance 
number  can  be  used  to  measure  the  goodness  of  the  residualization  and  for 
selecting  variables  to  be  residualized. 

The  two  options  of  the  reduction  algorithm  presented  above  are  implemented 
in  a  subroutine  called  REDUCE. 
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Shuffling  (Reordering  of  States) 


Denoting  the  reordered  (shuffled)  system  variables  by  x^,  u^  and  r^  the 
relation  between  reordered  variables  and  the  original  variables  is  expressed 

by 


x  =  P  x 

P  x 

u  =  P  U 

p  u 


(57) 


r  =  P  r 
P  r 


where  P  .P^and  Pr  are  the  shuffling  matrices  (obtained  internally  by 
software  from  the  shuffling  data  provided  by  the  user). 


Substituting  (57)  into  (48)  we  get 


x  =  A  x  +  B  u 

p  P  P  P  P 


r  =  C  x  +  Du 
p  p  s  p  p 


(58) 


where 

A  = 

P  AP_1 

P 

X 

X 

B  = 

P  BP_1 

P 

X 

u 

C  * 

P  cp"1 

P 

r 

X 

D  = 

P  DP  ‘ 

P 

r 

u 

(59) 
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Note  that  the  operations  described  by  (59)  to  obtain  A  ,  B  ,  C  and  D 

P  P  P  p 

are  merely  row  and  column  operations  (i.e. ,  permutations)  on  A,  B,  C 
and  D. 


The  shuffling  algorithm  is  implemented  in  the  subroutines  SHUFF,  SHUF  1 
and  SHUF  2  (SHUF  2  shuffles  the  name  list  table). 


Scaling  of  States 


Denoting  the  rescaled  variables  by  x  ,  u  and  r  ,  we  have 

S  S  S 


x  =  S  x 
s  x 


u  =  S  u 
s  u 


r  =  S  r 
s  r 


where  Sx,  S^,  and  Sr  are  the  scaling  matrices  (obtained  internally  by 
software  from  the  scaling  data  provided  ty  the  user). 


As  before,  the  scaled  system  is  described  by 


x  =  A  x  +  B  u 
s  s  s  s  s 


r  =  C  x  +  D  u 
S  S  S  s  s 


where 


SECTION  IV 


SYSTEM  PERFORMANCE  MODELING 


In  this  section  we  briefly  review  specific  performance  modeling  procedures 
for  the  ALDCS  design  process.  First,  general  performance  measures  are 
stated  for  completeness.  Subsequently,  the  steady-state  modeling  procedure 
is  described  in  detail.  This  is  followed  by  the  frequency  domain  modeling 
in  the  form  of  a  system  matrix.  Finally,  the  closed  loop  response  modeling 
is  given  with  output  feedback. 

GENERAL  PERFORMANCE  MEASURES 

As  illustrated  in  Section  II,  general  performance  measures  for  design 
and  analysis  can  be  listed  as: 

•  Poles  and  zeros 

•  Frequency  response  and  phase  and  gain  margins 

•  rms  response  to  turbulence  and  random  pilot  inputs 

•  Power  spectral  density  and  power  in  a  given  band 

•  Time  response 

Detailed  models  for  the  development  of  these  performance  measures  are 
given  in  References  1,  31,  and  50,  In  addition  to  these,  ALDCS  design 
requires  the  steady-state  values  of  time  response  with  step  inputs  for 
prescribed  steady-state  specification.  In  the  following  subsection,  we 
give  an  approach  for  this  requirement. 


STEADY-STATE  RESPONSE  MODELING 


Consider  a  system  described  by 

x  =  Ax  +  Bu. 

i 

r.  =  Cx  +  Du. 
x  1 

where  x,  r.,  and  m  are  the  state,  response,  and  input  vectors  of  sizes 
nx,  nr,  and  nu,  respectively. 

As  previously  discussed,  the  maneuver  load  control  calculations  involve 
computation  of  steady-state  control  surface  deflections  for  prescribed 
steady-state  responses.  For  the  steady-state  response  modeling,  the 
input  vector  m  in  general  consists  of  three  parts 

u.  =  col  [u  ,  u  ,,  u} 
i  o  d 

where 

u  =  set  of  inputs  with  prescribed  steady-state  values, 

ud  =  set  of  inputs  defined  by  prescribed  interconnections,  and 

u  =  set  of  inputs  with  unknown  steady-state  values. 

The  response  vector  is  also  divided  into  two  parts,  namely 

rg  =  a  set  of  specified  outputs  with  prescribed  steady-state  values  r^ 
r  =  total  outputs 


(63) 

(64) 
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The  input  vector  u.  and  the  specified  response  vector  r  can  be  represented 

by  the  steady-state  interconnection  quadruple  (P  ,  Q  ,  R  ,  S  )  as  follows: 

s  s  s  s 


u.  =  P 

1  s 

r.  +  Q  u 
x  s 

(65) 

r  =  R 

r.  +  S  u 

(66) 

s  sis 

Figure  14  shows  the  input/output  definition  for  the  steady-state  response 
calculations. 


Figure  14.  Interconnection  Model  for  Specified  Steady-State 
Inputs  and  Outputs 


The  set  of  equations  defining  the  steady-state  system  is  given  oy: 


Ax  +  Bu.  =  0 
i 

r.  =  Cx  +  Du. 
l  i 

u.  =  P  r.  +  Q  u 
1  sis 

r  =  R  r.  +  S  u 
s  sis 


r  = 


r 

o 


(67) 

(68) 

(69) 

(70) 

(71) 


r. 

i 
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•with  the  unknowns  x,  and  u.  When  r.  and  u.  are  eliminated  from  the  above 

1  l 

set  as  described  in  Reference  31,  we  obtain  the  following  set  of  equations 
in  terms  of  x  and  u. 


(72) 

Ax  +  Bu  =  0 

Cx  Du  =  r 

(73) 

o 

A  A 

r  =  Cx  +  Du 

(74) 

The  solutions  of  (72)  and  (73)  yield  the  required  steady-state  input  for  the 
prescribed  steady-state  output.  Equation  (74)  provides  the  steady -state 
value  of  the  total  response  vector. 

In  case  no  response  steady-state  value  is  specified  and  all  steady-state 
inputs  are  given.  Equation  (67)  is  directly  solved  for  x  and  the  total 
response  is  computed  from  (66). 

As  an  example,  consider  ALDCS  design  for  the  C-5A.  For  the  stick  per 
G  specification,  the  outboard  elevator  deflection  is  to  be  computed  for  1G 
maneuver  with  SAS  system  and  no  aileron.  The  steady-state  interconnection 
diagram  for  this  case  is  given  in  Figure  15. 


Figure  15.  Steady-State  Interconnection  Diagram 


The  input  and  specified  response  equations  are  given  by: 


6  =  0 
a 


ei 


KSAS  q  + 


eo 


6  =  u 

eo 

qUo 

acg  g 
where 

q  =  pitch  rate 

u  =  forward  velocity 

o 

g  =  acceleration  due  to  gravity 


(75) 


Using  these  equations,  the  interconnection  quadruple  is  formed  and  combined 
system  quadruple  (A  B  C  0)  is  obtained.  Finally  (77)  and  (78)  are  solved 
for  the  required  u  =  6  given  the  steady-state  value  a  =  -1G.  (See  Refer- 
ence  51,  Figure  18,  for  definitions  of  variables.) 


The  Honeywell  Software  which  implements  this  algorithm  is  called  Program 
SSK.  This  program  is  not  integrated  into  KONPACT  System  due  to  present 
scope  of  this  contract. 

FREQUENCY  DOMAIN  MODEL  FOR  OVERALL  SYSTEM 

(76) 

(77) 


Consider  the  system  described  by 
x  =  Ax  +  Bu 
r  =  Cx  +  Du 
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where  x,  r,  and  u  are  the  state,  response,  and  input  vectors  of  sizes 
NX,  NR,  and  NU,  respectively.  The  s-plane  (or  frequency  domain) 
description  of  the  system  in  the  FLEXSTAB/LSA  form  is  given  by 


~x  " 

[c(si] 

r 

u  _ 

where  C(s)  is  called  the  system  matrix  (Reference  2).  For  systems 
described  by  quadruple  (A,B,  C,  D),  it  takes  the  form  of 

C(s)  =  CjS  +  Cq 


The  steps  to  compute  Cj  and  Cq  from  the  quadruple  data  (A,  B,  C,  D) 
are  given  below: 


Taking  the  s-transform  of  Equations  (76)  and (77)  we  get: 

(si  -  A)  x  -  Bu  =  0 
r  -  Cx  -  Du  =  0 

Equations  (80)  and  (81)  can  be  combined  and  written  as  follows: 

si  -  A  !  0  1  -B 

1  I 

- t - j  -  " 

-C  'll  -D_ 
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(80) 


(82) 


***** 


Hence  and  are  given  by: 


1 

I  1 

0 

1 

!  0 

-  -4 

—  - 

T  ~  ~ 

0  ! 

0 

1  0 

-A  1 
1 

1  0 

1 

1  -B 

1 

1 

-C  ! 

1 

1  1 

1 

1  -D 

(£ 


( 


CLOSED-LOOP  MODEL  WITH  OUTPUT  FEEDBACK 


Consider  the  design  model  of  the  system  described  by 
x  =  Ax  +  Bu 

r  =  Cx  +  Du 


where  x,  r,  and  u  are  the  state,  response,  and  input  vectors  of  sizes 
NX,  NR,  and  NV,  respectively. 

The  input  and  output  vectors  are  in  the  form  of 
u  =  col  {uc,  T)3 

r  =  col  {r  ,,  r  } 
d  m 

where 


u  =  control  input  vector 


T|  =  disturbance  input  vector 


r^  =  design  response  vector 

r^  -  measurement  response  vector 

Equations  (85)  and  (86)  can  be  written  in  terms  of  these  subvectors  as 
follows 


• 

X 

=  Ax  +  B  u  +  B„  T\ 

1C  Cj 

(89) 

rd 

-  CjX  f  DjUc 

(90) 

r 

m 

*  C2X 

(91) 

The  closed-loop  equations  are  obtained  by  defining  the  output  feedback 
control  as: 


u  =  u  +  Kr 
c  m 

Substitution  of  (92)  into  (89) and  (90)  yields 
x  =  (A  +  Bj  KC2)  x  +  BjU  +  B2  7) 
rd  =  (Cj  +  D1  KC2)  x  +  Dj  u 


(92) 


(93) 

(94) 

(95) 


The  corresponding  closed-loop  quadruple  is  given  by 


A  =  (A  +  B  K  C_) 
1  ^ 

B  =  B  =  (B1  B2) 
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(96) 

(97) 


SECTION  V 


ACTIVE  CONTROL  SYNTHESIS  PROCEDURE 

INTRODUCTION 

This  section  presents  a  design  procedure  for  the  longitudinal  Active  Load 
Distribution  Control  System  (ALDCS).  First,  applicable  performance 
objectives  are  described  in  general.  Then  specific  ALDCS  design  goals 
are  stated.  Subsequently,  a  controller  configuration  for  the  Active  Load 
Distribution  Control  System  is  given.  This  is  followed  by  a  brief  descrip¬ 
tion  of  full  state  optimal  controller  design  and  simplified  optimal  controller 
design.  Finally,  the  ALDCS  design  results  us’ng  the  FLEXSTAB  C-5A 
residual  elastic  math  model  as  documented  in  References  5  and  G  are  presented. 

PERFORMANCE  OBJECTIVES  FOR  ACTIVE  CONTROL  SYNTHESIS 

The  goal  of  a  Control  Configured  Vehicle  (CCV)  design  is  to  improve  the 
performance  of  an  aircraft  using  active  control.  The  CCV  concepts  include 
the  following  areas  (Reference  13). 

•  Improved  handling  qualities 

•  Flight  envelope  limiting 

•  Reduced  static  stability 

•  Gust  acceleration  reduction 

•  Maneuver  load  control 

•  Active  control  of  structural  modes 
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In  the  following,  each  area  is  discussed  briefly  for  the  longitudinal  active 
control  design.  For  detailed  lateral  specifications,  see  Reference  13.  The 
enforcement  of  each  criteria  is  achieved  by  properly  weighted  response 
vectors  as  described  later. 

Improved  Handling  Qualities  (Handling  Quality  Control) 

The  longitudinal  axis-handling  qualities  are  specified  in  various  ways. 

The  major  ones  are  C*  model  response,  the  short  period  pole  locations, 
and  the  Time  Response  Parameter  (TRP). 

The  lateral  directional  handling  qualities  include  good  roll  rate  response, 
improved  Dutch  roll  damping,  and  good  turn  coordination  over  a  wide 
variation  in  angle-of-attack. 

Time  Response  Parameter  (TRP)  Criteria — The  command  performance 
is  specified  for  a  step  stick  input  in  terms  of  normal  acceleration  at  the 
c.g.  and  pitch  rate.  A  figure  of  merit  called  Time  Response  Parameter 
(TRP)  (see  Reference  15)  is  defined  to  measure  the  command  response 
(Figure  16).  The  requirement  is  for  a  TRP  below  a  specified  number. 


This  criterion  is  generally  applicable  to  any  system  order  and  degree  of 
linearity.  It  appears  to  correlate  reasonably  well  with  existing  specifica¬ 
tions,  e.g.,  MIL-F-8785B(ASG),  without  many  of  the  interpretation 
problems  of  the  latter.  It  is  amenable  to  both  computerized  performance 
evaluation  and  pilot  response  assessment.  Furthermore,  it  deals  directly 
with  input/output  relationships  which  are  significant  to  the  aircraft  function 
as  a  weapon  system. 
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Figure  16.  Definition  of  Time  Response  Parameter 


There  are  certain  difficulties  with  the  TRP  in  its  current  form,  however 
which  merit  consideration.  (See  Reference  14.) 


C*  Criterion — The  C*  criterion  is  an  example  of  specifying  short  period 
handling  qualities  in  terms  of  aircraft  parameters  familiar  to  a  pilot. 

The  concept  implicitly  includes  the  traditional  short  period  frequency  and 
damping  requirements  but  is  more  general  in  its  application.  The  usual 
definition  for  C*  is 


C*  =  K  N_  +  1^0  +  K  9 


a  z 


where  Ka>  K^»  Kc  are  dimensional  constants.  The  0  term  represents 
the  normal  acceleration  increment  at  the  pilot's  location  caused  by  the 


y 

I 


moment  arm  from  vehicle  center-of-gravity.  Therefore,  the  expression 
can  be  written  as 

U 

C*  =  N  +  K  =  — 

Z  Pilot  CO  * 

where  Kcq  =  "crossover"  velocity.  The  steady-state  perturbation  relation 

between  q  and  N  is 
z 

N 

z 

qgs  =  *  UQ  =  forward  velocity 

o 

The  velocity  at  which  the  contribution  of  pitch  rate  equals  the  contribution 
of  Nz  to  the  C*  response  is  the  so-called  crossover  velocity. 

The  C*  criterion  for  flight  control  has  evolved  because  it  allows  the 
designers  to  control  one  response  with  one  forcing  function  (the  elevator). 

At  high  dynamic  pressures  the  elevator  produces  primarily  normal  accelera¬ 
tion,  and  at  low  dynamic  pressures  it  produces  a  composite  variable  that 
is  significant  at  all  flight  conditions.  (References  52,  53) 

The  handling  qualities  can  be  summarized  as: 

•  The  dominant  short  period  frequency  as  excited  by  a  sharp-edged 
gust  shall  have  a  minimum  damping  ratio  of  0.  3 

•  For  a  step  pilot  input,  the  time  response  shall  meet  the  C* 
envelope  of  Figure  17  where  the  categories  are  defined  as: 
(Reference  13) 
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I  Optimum  response  (aerial  combat,  etc.) 

II  Not  as  critical  (refueling,  cruise,  etc. ) 

III  Categories  for  conditions  not  covered  by  1,  2,  4 

IV  Power  approach 

For  details  see  References  52  and  53. 


Figure  17.  Typical  C*  Envelope 

•  The  pitch  Command  Augmentation  System  (CAS)  will  produce 
a  steady-state  short  period  control  stick  gradient  consistent  with 
MIL-F-8785  requirements  (Reference  16).  A  C*  feedback 
permits  the  stick  gradient  requirements^©  be  met  without  scheduling 
the  stick  gain  (Kf).  This  is  illustrated  in  Figure  18  for  several 
points  representing  extremes  in  the  F-8C  data  (note  -  *>  -z  ) 
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n/a  (q/rad) 


Figure  18.  Typical  Pitch  Stick  Gradient 


Short  Period  Pole  Locations --In  this  criteria  the  natural  frequency  (ou^) 
and  the  damping  factor  (£)  of  the  short  period  dynamics  are  specified. 

It  is  the  simplest  criteria  for  handling  quality  specification. 

Flight  Envelope  Limiting  (Boundary  Control) 

The  flight  envelope  limiting  controller  design  involves  two  steps.  First, 
control  laws  are  defined  for  limiting  certain  variables  of  the  aircraft 
(i.  e.,  or)  during  maneuvers.  Second,  a  method  of  transition  (i.  e. ,  mode 
switching)  between  boundary  control  laws  and  normal  control  augmentation 
system  is  developed  (Reference  13). 

R  educed  Static  Stability  (RSS  Control) 


For  a  conventionally  designed  aircraft,  static  stability  and  acceptable 
handling  characteristics  are  obtained  through  aerodynamic  design  which 
includes  proper  location  of  the  c.g.  This  is  shown  in  Figure  19. 


t  Wing  and  Fuselage  Lift 


Figure  19.  Conventional  Static  Stability 
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In  maneuvering  subsonic  flight  and  in  supersonic  flight  this  usually  results 
in  significant  tail  down  loads  to  provide  the  required  moment  balance  for 
the  aircraft. 

If  a  high  authority  feedback  control  system  is  used  to  provide  artificial 
stability,  then  the  unaugmented  aircraft's  longitudinal  static  stability  can 
be  relaxed.  This  enhances  the  maneuvering  capability  of  the  aircraft 
by  reducing  the  drag.  Figure  20  shows  this  case. 


Figure  20.  Relaxed  Static  Stability  A/C  with  Active  Control 


The  reduced  static  stability  controller  (M^  augmenter)  restores  the  stability 
lost  in  the  CCV  aircraft  due  to  shifting  the  c.g.  aft.  It  basically  consists 
of  additional  pitch  rate  feedback  to  the  elevator. 


Gust  Acceleration  Reduction 


Reduction  of  the  aircraft  accelerations  due  to  wind  gusts  without  deteriorating 
the  response  to  pilot  commands  normally  enhances  the  stability  of  the  aircraft 
and  results  in  improved  mission  performance.  It  also  improves  the  ride 
qualities  for  the  pilot  as  well  as  enhancing  his  ability  to  perform  precision 
tracking  tasks.  It  is  well-known  that  the  conventional  elevator  together  with 
direct  lift  force  producers  (canards,  flaps  and  symmetric  ailerons)  can  be 
effective  in  reducing  gust  induced  accelerations.  Combining  these  control 
surfaces  with  the  elevator  controller  to  produce  direct  lift  for  gust  alleviation 
is  an  important  mode  in  active  control  laws. 

Maneuver  Load  Control  (Steady-State  Load  Relief  Control) 

For  transport  A/C  the  reduction  of  the  wing  root  bending  moments  during 
maneuvering  flight  (i.e.,  1  incremental  g  pitch  up)  is  described  to  alleviate 
structural  load  and  fatigue  of  the  wing.  This  reduction  in  wing  bending  is 
obtained  by  shifting  the  wing  lift  distribution  center  of  pressure  inboard  as 
shown  in  Figure  21. 


Figure  21.  Transport  Aircraft  Ideal  Lift  Distribution 


A  Maneuver  Load  Control  System  (MLC)  positions  control  surfaces  with 
steady-state  accelerations  to  redistribute  the  loading  on  a  wing.  For  fighter 
type  A/C,  the  objective  is  to  redistribute  the  wing  loading  to  reduce  drag 
during  high-g  maneuvers  (Reference  13). 

Structural  Mode  Damping  (Flexure  Control) 

Active  control  is  used  to  regulate  the  response  of  flexure  modes  to  turbulent 
air  or  pilot  commands  for  reducing  fatigue  damage.  These  can  be  extended 
to  damping  of  flexure  modes  for  flutter  suppression  also. 

By  placing  the  sensor  and  control  force  producers  near  one  another  the 
sensed  aeroelastic  forces  can  effectively  be  cancelled  out  by  the  control 
surfaces  to  lower  local  acceleration  (References  18  through  21). 

ALDCS  DESIGN  GOALS  -  FORMULATIONS  AND  PROCEDURES 

Active  Lift  Distribution  Control  System  (ALDCS)  Design  Goals 

One  of  the  objectives  of  this  contract  was  to  repeat  Honeywell's  ALDCS 
design  (Reference  4)  using  the  FLEXSTAB  generated  C-5A  A/C  data  and 
KONPACT  software.  The  cruise  flight  condition  was  selected  for  this 
demonstration  as  shown  in  Table  2. 

Full  State  Quadratic  Design  Formulation 

The  theory  and  numerical  techniques  used  in  quadratic  design  are  well 
documented  in  other  sources  (References  12,  23,  24,  and  25)  and  therefore 
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will  not  be  repeated  here.  We  shall  be  concerned  here  primarily  with 
formulation  of  the  problem.  The  system  description  has  the  form 


x  =  Ax  +  B^u  +  B2T| 


r^  =  CjX  +  DjU 


r2  =  C2X 


(100) 

(101) 

(102) 


where 


x  =  state  vector  (including  rigid-body  states,  actuator  and  servo  states, 
flexure-mode  states,  sensor  states,  model-following  states,  and 
wind  states) 

u  =  control  input  vector 

11  =  unit-variance  white  noise  vector 

Tj  =  design  response  vector 

r9  =  measurement  vector  (consists  of  actual  and  complementary 

Ci 

measurements ) 


We  wish  to  find  a  time-invariant  controller  of  the  form 


u  =  K*r2  =  K*C2x 


(103) 


which  minimizes  the  performance  index 


J  =  E{r^  Q  r^  =  tr[Q  Rj] 


(104) 
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where  Q  is  a  symmetric  weighting  matrix,  and  is  the  response  covariance 
matrix  given  by 

R1  =  (C1  +  DiK  '<C2)  X  (C1  +  DiK>"C2)/  (105) 

and  X  is  state  covariance  matrix  given  by  the  solution  to  Equation  (108). 

Figure  22  shows  the  block  diagram  of  the  overall  system. 


(Design  Response) 


(Actual  and 
ComDlementarv 
Measurements) 


Figure  22.  Optimal  System  Block  Diagram 


The  optimal  gain  on  the  full  measurements  is  given  by 

K*  =  -(DJ  QDj)'1  (D^  QC1  +  P)  X  (C2  X  C^)"1  (106) 
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where  P  (Riccati  matrix)  and  X  (covariance  matrix)  are  determined  from 


(A  +  Bx  K*C2)'  P  +  P(A  +  BjK^)  +  (Cj  +  DjKC^'  Q(C  +  D  K*C  >  =  0 

(107) 

(A  +  B^Cg)  X  +  X(A  +  B1K*C2)+  B2  =  0  (108) 

When  is  invertible  (for  full  measurement  feedback),  the  solution  of 
Equation  (106) 

K*  =  -(DJ  Q  Dj)"1  (D’  QCj+BJ  Pjc"1  (109) 

and  does  not  depend  on  the  covariance  X.  For  limited-measurement  feed¬ 
back,  Cg  is  not  invertible  since  the  number  of  measurements  is  generally 
fewer  than  the  number  of  states,  so  that  K*  will  be  a  function  of  both  P  and  X. 

Simplified  Quadratic  Design  Formulation 

The  design  procedure  described  in  References  12,  23,  and  25  is  used 
to  simplify  the  full  state  controller.  In  this  procedure,  the  measurement 
gains  are  written  as  a  function  of  a  scalar  parameter,  X,  such  that 

u  =  K*(X)  r2  (110) 

where 

K*(X)  =  Kj(X)  +  X  K2  0  £  X  £  1 

The  starting  point  (X=l)  is  fou:  d  by  using  the  optimal  state  feedback  gains 
and  the  measurement  matrix  (corresponding  to  actual  measurements 


•--.-  .•ft*'i«f!»-V’w-  •»**«  •' 


augmented  with  complementary  measurements  so  that  exists) 

K*(l)  =  KC”1  (111) 

The  measurement  constraints  are  applied  gradually  by  stepping  X  to  zero, 
thereby  reducing  gains  on  the  complementary  measurements.  The  matrix 
K^(0)  is  the  fixed-form  solution  and  has  the  gain  structure  desired. 

This  procedure  of  "backing  off"  from  the  state  feedback  controller  is 
illustrated  in  Figure  23. 


Figure  23.  Optimal  Control  with  Measurement  Constraint 


The  same  quadratic  performance  index  as  of  the  full  state  feedback  is 
minimized  using  the  gains  after  each  stepping  down  of  X  parameter. 
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C°nstrained  Quadratic  Design  Formula tion 


As  will  be  shown  later,  the  maneuver  load  control  requirements  generate 
an  equation  in  the  form  of 


where 


tate  of  the  overall  system  at  one  incremental  g 


required  steady-state  aileron  deflection 
by  specified  amount 

aileron  feedback  gains 


to  reduce  bending  moment 


Figure  24.  Optimal  Control  with  Gain  Constraints 


This  problem  can  be  handled  the  same  way  by  performing  a  transformation 

A 

in  the  gain  space  and  stepping  down  of  \  parameter  on  and  optimizing 
with  K^. 

Response  Selection  Procedure 

Quadratic  design  formulation  shows  that  design  problem  is  essentially 
reduced  to  defining  proper  responses  (response  selection)  for  the  cost 
function  and  finding  proper  weights  (weights  selection)  for  the  performance 
achievement. 


These  two  steps  require  design  experience  and  insight  to  the  physics  of  the 
problem  as  well  as  familiarity  with  the  synthesis  software. 

Response  selection  for  design  depends  on  the  design  criteria  to  be  enforced. 
Converting  a  set  of  design  criteria  into  a  minimization  of  a  set  of  responses 
is  one  of  the  fundamental  steps  in  optimal  control  synthesis. 

In  the  following,  we  will  discuss  typical  cases  briefly  to  demonstrate  the 
mechanics  of  the  response  selection. 


Response  Selection  for  Enforcing  Handling  Quality  Criteria--The  C-5A 
aircraft  is  augmented  with  a  simple  Stability  Augmentation  System  (SAS) 
to  enhance  the  handling  qualities  of  the  aircraft.  It  consists  of  pitch  rate 
feedback  to  the  inboard  elevator  input  in  the  form  of 


6  . 
eic 


KSASq 


(113) 
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One  way  of  incorporating  the  effect  of  this  SAS  into  optimal  handling  quality 
controller  design  with  an  aileron  input  is  to  define  the  inboard  elevator 
residual  response  (deviation  from  ideal)  using  (113). 

r  ‘  •  6eic>  <114) 

Thus  by  varying  the  weights  on  r,  the  short  period  roots  can  be  indirectly 
controlled. 


Implicit  Model  Following 

When  handling  quality  criteria  specifies  the  desired  short  period  locations, 
as  in  C-5A  case,  the  model  following  error  rate  response  may  be  generated 
to  enforce  the  specified  pole  locations  in  the  following  way. 


The  rigid  body  (KB)  and  model  M  equations  are  given  by 


XRB  ~  ARB  XRB  +B6e 


XM  *  AM  XM  +  B6e 


The  rate  of  residual  is  defined  as 


r  =  dt  ^xrb  "  XM^ 


and  given  by 

r  =  ARB  xrb  ”  am  xm 


Assuming  that  rigid  body  and  model  states  are  close  together 


(115) 

(116) 


(117) 


(118) 


XM  ~  XRB 


(119) 


yields 


r  =  (ARB  "  am)  xrb 


(120) 


Let  ouM  and  be  the  natural  frequency  and  damping  factor  of  the  model 
short  period  modes.  The  corresponding  characteristic  equations  are 
given  as 


1,2  +  (2Cm  V  s  +  “m  °  0 


(121) 


In  terms  of  the  model  transition  matrix 


s  +  (tr  A.J  s  +  det  A  =  0 
M  M 


(122) 


Equations  (121)  and  (122;  provide  two  equations  to  determine  the  elements 
of  model  transition  matrix  A.^.  If  we  assume  that  the  w  equations  for  model 
and  vehicle  are  identical,  then  the  remaining  elements  of  A^  can  be  found 
from  (121)  and  (122).  Equation  (120)  is  the  response  equation  to  force  the 
desired  pole  locations,  with  the  coefficient  matrix  as  computed  above. 


If  symmetric  ailerons  are  used  together  with  the  inboard  elevators  for 
Manevuer  Load  Control  (MLC)  and  Gust  Load  Alleviation  (GLA),  then 
handling  quality  can  be  enforced  in  the  following  way. 

Figure  25  shows  the  block  diagram  of  the  rigid  body  part  of  the  flexible 
vehicle  dynamics  and  corresponding  handling  quality  model  which  generates 
existing  SAS  pitch  responses  (w^,  q^). 
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R.B.  Vehicle 
Dynamics 


"  1 


eiri^l  R.B.  Vehicle 


Dynamics 


Model 


Figure  25.  Response  Generation  to  Enforce  Handling  Quality  Criteria 


The  rigid  body  part  of  the  flexible  vehicle  equations  are  in  the  form  of: 


XRB  =  A11  XRB  +  A12  XBM  +  A13  XA  +  A14  XG 


where 


=  Rigid  body  states 


=  Bending  mode  states 


=  Control  surface  states 


=  Gust  states 


Pilot  states 


Handling  quality  open-loop  model  is  in  the  same  form  of  rigid  body  equations 
and  given  by 


XM  ~  A11  XM  +  A12  XBM  +  A13  XAM  +  A14  XG 


(124) 


where 


=  model  states 


XAM  ~  moc^el  control  surface  states 


Here  the  model  control  surface  states  are  specified  as 


6  . 
eiM 


KSAS  qM  +  Xp 


(125) 


This  assumes  that  the  SAS  aircraft  is  the  ideal  model.  The  residual 


response  rate  on  the  rigid  body  state  is  defined  as 


r  =  x  —  x 
RB  RB  M 


(126) 


Substituting  (123)  and  (124)  into  (126)  yields 


rRB  ~  A11(XRB  “  XM)  +  A 13  (  XA  ' 


(127) 


Assuming  that 


XRB  *  XM 


(128) 
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Equation  (127)  yields  the  residual  response  rates  in  terms  of  the  linear 
combination  of  states  and  the  system  coefficients 
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(129) 


Thus  by  varying  the  weights  on  r..^,  the  handling  quality  is  enforced  in 

K  ri 

optimal  ALDCS  design. 

Explicit  Model  Following 

In  the  preceding  analysis,  the  response  equations  contain  only  the  original 
system  states.  The  model  states  are  discarded  using  the  rates  of  the  residuals 
and  the  assumption  that  model  and  system  states  are  approximately  the  same. 
This  technique  is  known  as  the  implicit  model  following.  It  has  the  benefit 
of  not  increasing  the  system  order  and  not  yielding  high  gains  in  the  controller 
design. 

In  the  following,  we  present  for  completeness  another  form  of  response 
development  to  enforce  handling  qualities  which  is  known  as  explicit  model 
following. 

Figure  26  shows  handling  quality  controller  which  produced  improved  short 
period  response  with  an  explicit  model  (Reference  13).  The  integrator 
in  the  controller  enforces  the  following  equation  in  the  steady-state 

‘ =  yM  -  y  = 0 
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cos: 

3 

+ 

tn 

s: 

cm  s: 

Is: 

3 

kj' 

<NJ 

4* 

CM 

tn 

■< 


where 


v  =  Model  response  to  command 
J  M 

y  =  Blended  pitch  rate  and  lagged  normal  acceleration 

The  design  response  to  enforce  the  desired  handling  quality  is  in  the  form 
of : 

• 

r  =  col  (e,  x  ,  6  ) 
e  e 

which  contains  the  model  following  error  and  its  integral.  It  also  contains 
actuator  rate  response  to  limit  the  bandwidth  of  the  actuator  due  to  feedback. 

Figure  27  shows  the  general  structure  of  controller  with  explicit  model 
(Reference  23).  The  actual  model  in  the  control  law  to  generate  model 
response  x^  increases  the  dynamical  order  of  the  controller  and  generally 
yields  higher  gain  values  as  compared  to  controller  with  implicit  model. 

( - 1 

Model -Following  i 

Controller  I 


Figure  27.  Explicit  Model- Following  Structure 


The  control  law  is  in  the  form  of 

u  =  K*  +  K*..x  +  K*  x 

x  MM  P  p 

where 


x  =  plant  state  including  vehicle,  actuator,  sensor  and  gust  filters 


=  model  states 
=  pilot  input  states 


Response  Selection  to  Enforce  Maneuver  Load  Control — Maneuver  load 
control  involves  computing  bending  moment  response  of  free  A/C  for  1G 
normal  acceleration  pitching  maneuver  at  the  level  flight  and  determining 
control  surface  positions  to  reduce  the  bending  moment  as  specified  by  the 
specification  (see  Figure  28). 

The  normal  acceleration  is  given  by 

a  =  w  -  U  q  (132) 

eg  o  H 

where  Uq  is  the  cruise  speed  and  w  =  o  for  level  flight. 


Figure  28.  Maneuver  Load  Control  Surface  Position  Calculation 


First  pilot  stick  input  to  free  A/C  as  shown  in  Figure  28  is  calculated  to 
produce  a  steady-state  pitch  rate  corresponding  to  1G  normal  acceleration. 
The  value  of  the  speeificU-bending  moment  is  obtained  at  this  maneuver. 

Subsequently,  the  pilot  input  for  1G  maneuver  with  SAS'd  A/C  is  calculated 

to  enforce  the  stick/g  requirement.  This  determines  the  steady-state 

deflection  for  the  outboard  elevator  6  .  Now  the  aileron  and  inboard 

eo 

elevator  deflections  are  computed  to  obtain  1G  maneuver  and  reduced 
bending  moment  response,  with  the  fixed  value  of  6  . 

After  having  determined  the  steady-state  surface  positions,  the  next  task 
is  to  enforce  these  positions  by  properly  selected  controller  configuration. 


Maneuver  load  control  can  be  enforced  either  by  integral  control  or  by  a 
high-passed  aileron  control.  The  integral  control  enforces  the  following 
aileron  deflection  equation  at  steady-state. 


n 

a 


=  0 


Figure  29  shows  the  MLC  configuration  with  integral  control. 


(133) 


The  integral  control  has  one  obvious  drawback  -  it  may  require  a  sensor  to 
measure  aileron  position  (Reference  8).  High  passed  aileron  control,  as  shown 
in  Figure  30,  washes  out  all  steady-state  aileron  commands  except  for  MLC. 

CONTROLLER  CONFIGURATIONS 


Handling  quality  controller  with  explicit  model  and  state  feedback  corresponding 
to  Figure  26  is  shown  in  Figure  31.  The  desired  reduced  form  is  shown  in 


6  =  Aileron  actuator  input 

clC 

Kj  =  Gain  integral  of  error  state 

£  =  Error 

x^,  =  Error  state 

x  =  Total  system  state 
K  =  State  feedback  gain  matrix 


Figure  29.  Integral  Control  Configuration  to  Enforce  MLC 


Figure  30.  High-Passed  Aileron  Control  Configuration  to  Enforce  MLC 


Figure  32.  Boundary  controller  require  isolation  of  a  common  proportional 
plus  integral  function,  as  given  in  Reference  13.  Gain  scheduling  is  also 
facilitated  by  this  configuration.  The  factorization  process  is  not  unique. 
One  such  factorization  is  given  in  Reference  13. 

Figure  33  shows  the  configuration  used  in  Reference  4.  Design  problem  is 
to  determine  the  coefficients  of  the  second  order  maneuver  load  control 
and  gust  alleviation  filter  and  the  four  feedback  gains  shown  in  the  diagram. 
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Handling  Quality  Controller  Configuration  with  Explicit  Model 


Figure  32.  Desired  Form  of  Handling  Quality  Controller  Configuration 
(Proportional  and  Integral  Control  Factorized) 


This  configuration  provides  steady-state  aileron  deflection  proportional  to 
steady-state  normal  acceleration  as  given  by 

(134) 


6a  KSS  A21R 


The  Kgg  is  obtained  from  the  steady-state  maneuver  load  control  require¬ 
ments  as  described  previously.  In  addition,  this  filter  shapes  the  frequency 
response  of  the  aileron  loop  to  provide  gust  load  alleviation. 


The  decomposition  of  this  filter  into  a  lagged  acceleration  and  an  integrator 
is  shown  in  the  state  diagram  as  given  by  Figure  34.  Table  8  provides 
a  relationship  between  the  filter  coefficients  and  the  optimal  feedback  gains. 

Again  this  decomposition  is  not  unique,  and  can  be  accomplished  in  different 
ways.  For  instance,  the  filter  decomposed  as  shown  in  Figure  34  emorces 

the  relation  given  by  Equation  134  whereas  the  decomposition  given  in  Figure  38 
enforces  an  equivalent  relation  given  by 

6  6a  =  KAF  (134a) 
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Table  8.  Relation  of  Filter  Coefficients  to  Controller  Gains 
(as  implemented  in  Figure  34) 


In  both  configurations  additional  states  are  introduced  into  the  control  law. 
These  states  are  included  in  the  response  set  so  that  they  can  be  weighted 
for  proper  overall  response. 

The  lagged  normal  acceleration  (A21RL)  time  constant  is  chosen  to  be  the 

same  as  aileron  actuator  (6  )  time  constant.  This  selection  eliminates  the 

a 

need  for  the  aileron  actuator  sensor. 

Past  Honeywell  work  (Reference  4)  indicates  that  the  third  bending  mode 
damping  is  difficult  to  control  during  the  ALDCS  design.  To  facilitate 
design  and  reduce  damping  factor  sensitivity  to  gains,  one  may  include  a 
roll-off  filter  in  the  elevator  loop  as  shown  in  Figure  33. 

C-5A  ALDCS  EXAMPLE 

The  ALDCS  example  presented  here  is  based  on  Reference  4  and  utilized 
to  verify  the  KONPACT  software  results.  The  ALDCS  design  goals  are  shown 
in  Table  9.  The  C-5A  ALDCS  controller  design  via  the  optimal  control 


Table  9.  ALDCS  Design  Goals 


reverse  directions 


synthesis  technique  is  carried  out  in  three  steps:  a)  Design  Model  Genera¬ 
tion,  b)  Controller  Design,  and  c)  Performance  Evaluation. 


Design  Model  Generation 


The  detailed  steps  to  obtain  the  design  model  are  described  in  Section  VI 

of  Volume  III  (Reference  51).  Essentially  it  involves  converting  the  FLEXSTAB/ 

LSA,  15  mode,  residual  elastic  simulator  deck  data  into  state  space  form 

and  augmenting  it  with  the  gust  model  and  actuator  to  obtain  the  plant  model. 

These  are  done  by  the  software.  This  plant  model  is  subsequently  reduced 

by  residualization  procedure  (see  Section  VI)  and  the  reduced  model  is 

augmented  with  the  controller.  After  this,  design  responses  are  selected 

to  obtain  the  design  model.  Figure  35a  and  b  show  the  input/output  diagram 

of  the  design  model.  The  actuator  state  diagram  is  shown  in  Figure  36. 

The  form  of  the  gust  model  transfer  function  block  diagram  is  shown  in 

Figure  37.  Parameter  values  are  given  in  Reference  4.  The  ALDCS 

controller  state  diagram  used  in  Reference  4  is  shown  in  Figure  38.  In 

the  controller,  the  constants  and  KM1  are  determined  to  satisfy  the 

steady-state  maneuver  load  control  requirements.  The  steady-state 

response  computations  on  the  plant  model  are  illustrated  in  Table  10.  In 

the  first  case  we  set  6  =0  and  5  .  =  6  =  p.  The  value  of  p  is  determined 

a  ex  eo 

to  obtain  q  =  q  =  0.04377.  (For  consistency,  the  steady-state  values  are 

computed  at  +1  incremental  g,  i.e. ,  2g,  in  level  flight.  At  the  cruise  speed 

of  735  fps  this  corresponds  to  a  pitch  rate  of  0.04377  rad/sec.)  The 

corresponding  steady-state  value  of  is  also  obtained.  This  is  designated 

as  B..  In  the  second  case,  we  set  5  =  0,  5  .  =  p  +  0,5  q  and  6  =  p 

and  again  the  value  of  p  is  determined  to  obtain  q  =  q.  This  is  designated 

as  p.  In  the  third  case,  we  set  5  =  p  and  compute  6  and  6  .  to  obtain 

60  a.  61 

q  =  q  and  =  0.7  B^. 
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Design  Model  Generation  for  ALDCS  Controller  Design 
(C-5A  Cruise  Flight  Condition) 
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ACG/f 
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f  the  maneuver  load  control  gain  is  calculated  as 
=  6,0  S&/g  B  (6.0)  (-0.  2752) 
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(1/6.0)  (-1.0) 


Figure  38.  ALDCS  Controller  State  Diagram 
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The  steady-state  value  6  ^  is  used  lu  adjust  the  feedback  gains  of  the  ALDCS 
controller  in  the  end  to  meet  the  steady-state  requirements  for  maneuver 
load  control. 


Controller  Design 


The  ALDCS  controller  design  presented  here  is  repeat  design.  Figure  39 
represents  a  block  diagram  of  the  design  process.  The  design  response 
weights  are  shown  in  Table  11.  The  following  steps  describe  the  controller 
design. 


White  Noise 
Disturbance 
Inputs 


Control 

Inputs 


Overall  System 
Design  Model 


Design 

Response 


=n] 


Measurements 


Feedback 

Gains 


Figure  39.  ALDCS  Controller  Design  Process 
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Responses 

Weights 

Value 

MLCl 

0.  800E+01 

B1 

«2 

0. 100E-01 

T1 

,o 

co 

0. 100E-0P 

«4 

0 

B2 

0 

T2 

Q6 

0 

• 

6 

a 

Q7 

0.  500E+04 

B3 

0 

T3 

0 

• 

6  . 
ei 

yio 

0.600E+06 

B4 

Qli 

0 

T 

4 

Q12 

0 

6 

a 

Q13 

0 

B5 

Q14 

0 

T5 

Q15 

0 

6  . 
ei 

Q16 

0 

Bi 

Q!7 

0.  750E-13 

Q18 

0.  100E-10 

Table  11.  Design  Response  Weigh 
(C-5A  Cruise  Flight  Co 


Weights 


Value 


0.  100E-13 


0. 100E-11 


0.  200E-13 


0.  200E-11 


0.  800E-13 


0. 100E-10 


0.  200E-12 


0.  200E-10 


0. 100E+01 


0. 100E+01 


Full  State  Optimal  Controller  Design--Full  state  optimal  control  is  unique 
and  is  used  for  initial  studies.  It  determines  what  performance  can  be 
achieved  from  the  system  under  idealized  conditions. 

Using  the  overall  system  design  model  described  earlier,  the  optimal 
full  state  feedback  gains  are  computed  for  the  specified  quadratic  weights. 


Simplified  Optimal  Control  Design- -Starting  from  the  optimal  full  state 
feedback  gains,  simplified  feedback  gains  are  obtained  as  described  in 
References  25  and  4.  The  reduced  control  law  is  given  by  the  following 
equations  (see  Figure  38). 

u6  ■  <K1a>  6  a  +  (k1A21rl> A21RL  +  ,k1glaf>  glaf  (i35) 

a  a 


U6  =  (K2A21R>  A21R  +  <K2AFUS)  AFUS  +  ^TFUS1  TFUS  + 
ei 

(136) 

(K2p)P 

The  reduced  feedback  gains  are  shown  in  Table  12. 


Adjustment  of  Gain  K2p  and  Calculation  of  the  MLC  Filter  Coefficients -- 

The  gain  K2  obtained  in  the  previous  step  is  adjusted  so  that  the  steady - 
P 

state  maneuver  load  control  requirements  are  met  for  the  elevator  channel. 
From  the  steady-state  response  calculations 

6  .  =  -0.0226  rad  5  =  -0.0896  rad 

ei  eo 


Since  the  actuator  steady-state  gain  is  unity  (refer  to  Figure  38)  we  have: 
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Table  i2.  Reduced  Feedback  Gains  for  ALDCS  Controller  Design 
(C-5A  Cruise  Flight  Condition) 


Gains 

Values 

k1dela 

-7. 81200 

K1A21RL 

11.94000 

k1glaf 

1.96900 

K2A21R 

0.00256 

K2afus 

-0.06400 

K2tfus 

0.49040 

K2 

P 

-0.40150 

6ei  =  U«.  ’  V  +  P  ■  <K2A21R>  A21R  + 
ei  2 


(137) 


<K2afus’  afds  +  (K2tfus>  tfus  +  <K2*p  +  1-0)P 


Substituting  the  steady-state  values  for  1G  maneuver  and  solving  for  K2*^, 

we  obtain  K2*  =  0.178. 

P 


Table  13  shows  the  adjusted  reduced  gains.  Figure  40  shows  the  final  form 
of  the  ALDCS  controller  configuration  for  cruise  flight  condition.  The  final 
maneuver  load  control  and  gust  alleviation  filter  coefficients  are  obtained 
from  A21RL  and  GLAF  filters  (Figure  38)  by  setting  the  integrator  feedback 


i 


gain  AF  =  0.0.  These  calculations  are  shown  in  Table  14.  The  final  ALDCS 
controller  is  shown  in  Figure  41. 

Table  13.  Final  Gains  for  ALDCS  Controller  Design 
(C-5A  Cruise  Flight  Condition) 


Gains 

Values  5 

k1dela 

-7. 812000 

K1A21RL 

11.940009 

k1glaf 

1.969000 

K2A21R 

0.002565 

K9 

AFUS 

-0.064010 

K9 

TFUS 

0.490400 

K2  * 

P 

0. 178000 

*  K2  Gain  obtained  by  FFOC  is  subsequently  adjusted  to  satisfy  the 
steady-state  ALDCS  requirements  for  5 

ei 


Table  14.  MLC  and  GLAF  Filter  Coefficients 
(As  implemented  in  Figure  38) 


|  Filter 
Coefficient 

Equation 

Value 

Ai 

6(1-KW> 

52. 872 

A 

o 

36*k1glap 

70.  88 

Bi 

TCI 

A21RL 

11.  94 

B 

0 

KAF*K1 

aglaf 

19.  51 
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Maneuver  Load  Control  and  Aileron 

Gust  Alleviation  Filter  ALDCS  Off  Actuator 


Figure  40.  ALDCS  Controller  Configuration  for  Cruise  F.  C.  412301 


Aileron  Actuator 


ALDCS  Block  Diagram  for  Cruise  F.  C.  412301 


Performance  Evaluation 

The  detailed  performance  results  are  presented  in  the  next  section  Here 

the  variation  of  some  hey  performance  variables  as  w.  move  from 

state  feedback  to  simplified  feedback  (,e„  as  (l-X)goes  to  1,  is  presented. 

Figures  42,  43,  and  44  present  the  variation  of  quadratic  cost  and  the  sensor 
.  /,  ti  varies  from  0  to  1.  These  vanatxons  are  lair  y 

- ».  steP  « ,  > — 

large).  Note  that  the  aileron  feedback  gains  to  taking 

(1-k)  approaches  aero  and  then  suddenly  s  °°  “^  integrator  in  Figure  38. 

out  the  temporary  MLC1  integrator  and  putting  in 

Figures  45.  47.  49.  and  51  present  the  variation  of  rms  responses  <e 
g  .  ,  .  _  M  x\  varies  from  0  to  1.  And  Figure  46, 

e  .  Tl,  Bl)  due  to  wind  gust  as  (1  X)  v  #  . 

4g  50  and  52  present  the  variation  of  rms  responses  (eq,  ew>  Tl, 

,  ,  ■  from  0  to  1.  in  all  of  these  figures,  we  observe 

that  7uU '°tate  feedback  determines  what  performance  can  be  obtained  from 
the  system  under  ideal  conditions  and  the  simplified  feedback  is  obtainedby 

n  minimum  loss  of  the  performance  simplified 

note  that  the  final  ALDC.3  performance  is  slightly  different 
(reduced)  feedback  controller  performance.  (Final  ALDCS 

the  steady-state  retirements  completely  whereas  the  amplified  ee^iack 
controller  does  not.  This  is  because  of  the  small  feedback  around  he  0 
filter  (i.  e. ,  AF  f  0  in  Figure  38)  in  the  simplified  controller  as  well  the 
adjustment  of  pilot  gain  K2p  to  K*2p  for  meeting  steady-state  requiremen  . 

Figures  53,  54.  and  55  present  the  variation  of  the  closed-loop  roots  (rigid 
body.  firs,  flexure  mode,  and  third  flexure  mode)  as  (1-X)  vanes  from  1. 
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Figure  45.  Variation  of  RMS  Response  (For  Implicit  Model  Error 
Rate  for  q  (e  )  Due  to  Wind  Gust)  with  (l-X) 
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Figure  46 
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Figure  47.  Variation  of  RMS  Response  (For  Implicit  Model  Error 
Rate  for  w  (e  )  Due  to  Wind  Gust)  with  (1-X) 
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Figure  50.  Variation  of  RMS  Response  (For  Torsion  Moment  (Tl) 
Due  to  Pilot)  with  (1  -X) 
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SECTION  VI 


IMPACT  OF  MODELING  AND  MODEL  REDUCTION  PROCEDURES 
ON  CONTROL  SYSTEM  DESIGN  AND  PERFORMANCE 


Complete  representations  of  flexible  airplanes  requires  models  of  very 
high  order  (100  or  more).  These  complete  models  must  be  used  for  final 
analysis  and  verification.  Most  control  synthesis  can  be  determined  from 
much  lower  order  models.  The  subject  under  discussion  is  how  to  make 
the  best  low  order  approximation.  Best  cannot  be  well  defined.  Roughly, 
we  desire  to  synthesize  using  the  low  order  model  and  hope  that  only  minor 
tailoring  is  required  to  obtain  comparable  results  on  the  complete  model. 

In  the  following,  model  reduction  procedures  are  reviewed  briefly.  Then 
the  truncation  and  residualization  reduction  procedures  used  in  KONF'ACT-1 
are  described.  This  is  followed  by  a  detailed  comparison  of  Honeywell/ 
GELAC  and  FLEXSTAB  data  on  C-5A  aircraft  using  the  truncation  and 
residualization  operations, 

REVIEW  OF  MODEL  REDUCTION  PROCEDURES 

In  the  following,  closed  and  open-loop  approximations  are  discussed. 

Closed-Loop 


Considerable  progress  towards  developing  synthesis  and  analysis  of  the 
regulator  problem  for  high  order  systems  is  presented  in  References  32 
and  33.  The  plant  is  taken  to  be 


I 

H 

E 


■ 


Xx2  =  Vl  +A4X2  +  B2U 


(139) 


where  H  is  a  small  parameter.  If  X  =  0,  the  usual  residualization  approxima¬ 
tion  is  being  made. 


Synthesis  and  analysis  are  performed  after  invoking  the  quadratic  performance 
criterion 


J  =  J  {y'Qy  +  u’Ru}dt 
o 


(140) 


where 

y "  cixi +  C2X2 


R  is  positive  definite 
Q  is  positive  semi-definite. 


Theorems  and  examples  show  continuity  with  respect  to  X  as  X  -*  0  .  Refer¬ 
ences  32  and  33  present  the  only  results  found  on  the  closed-loop  approximation 
problem.  They  have  made  considerable  progress;  the  procedures  appear 
to  be  sound,  tractable,  and  capable  of  extension.  These  efforts  should  be 
continued  toward  resolving  whether  they  can  provide  a  better  solution  to 
control  synthesis  for  flexible  aircraft. 
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Open-Loop 

There  is  a  surfeit  of  open-loop  approximation  schemes  to  reduce  the  order 
of  large  systems.  Unfortunately,  they  have  not  been  ranked  (even  for  a 
simple  application).  Promising  schemes  for  applicability  to  the  flexible 
vehicle  are  discussed  below. 


Truncation  and  Residualization--With  truncation,  x^  and  x2  are  set  equal  to 
zero  in  Equation  (138);  Equation  (139)  is  omitted.  This  approximate  model 
for  high  order  systems  has  as  its  main  virtue  simplicity  in  application. 

For  a  given  reduced  order  model  it  is  usually  the  worst. 


With  Schwendler  and  MacNeal's  residualization  (Reference  46),  \  is  set 
equal  to  zero  in  Equation  (139).  This  approximation  yields  exact  values 
with  respect  to  the  steady-state  response  of  the  large  system.  Further 
details  on  these  are  given  in  Section  III  on  Model  Conditioning  and  in  this 
section. 


Stable  Partial-Pade  Moment  Matching--References  34  through  38  present  rather 
simple  computational  methods  for  approximating  the  system  H(s)  by  an 

a 

approximate  models  H^(s)  where 


and 


X(s)  =  H(s)  U(s) 

bs11  *  + - - - +b 

H<s)  '  ^ - - 

n  n  —  l 

as  +a„s  + - +a 

o  1  n 

A  k-1  A 

-  b  s  1  + . +b 

H„<s>  =  ~~g~ r-^i - - 


(141) 


(142) 


(143) 


as  +  a,s 


where 


X(s)  =  p-vector 

U(s)  =  r-vector 

b.  =  p  x  r  matrices 
1 

a.,  a.  =  scalars 
i  3 

k  <  n 

The  approximation  H^(s)  is  stable  if  H(s)  is  stable. 

Let  the  MacLaurin  expansion  of  H(s)  and  H^(s)  be 

H(s)  =  “0  +  o^s  +  a^s2  +  -  -  -  -  (144) 

H.  (s)  =  &  +  a  s  +  a_s2  +  - -  (145) 

k  o  1  2 

The  criteria  for  obtaining  reduced  system  is  to  match  the  coefficients  in 
(144)  and  (145)  starting  with  the  lower  degree  terms. 


Another  interpretation  of  (144)  or  (145)  is  available.  For  a  single-input 
single-output  system 

00 

a.(il)  (-1)1  =  JQ  tXh(t)  dt  (146) 

-  the  i**1  time-moment  of  h(t) 
where  h(t)  is  the  impulse  response  corresponding  to  H(s). 


It  is  seen  that  this  is  a  generalization  of  residualization:  exact  not  only  in 
the  steady-state  but  also  a  good  approximation  at  low  frequencies. 


Iterative  Matrix  Matching  By  Identification — Reference  39  presents  an 
identification  scheme  that  has  been  used  at  Honeywell  to  construct  lower 
order  approximations  to  high  order  models. 


x  =  Ax  +  Bu 

(147) 

y  =  Cx  +  Du 

(148) 

*A  A  A  A 

x  =  Ax  +  Bu 

(149) 

*"h 

Equation  (147)  is  the  known  plant  of  n“  order  with  prescribed  response 

th 

relationship  (148).  Equation  (149)  is  an  m  order  approximation  with  m<n. 

A  A 

Both  systems  are  forced  by  a  prescribed  u.  A  and  B  are  determined  to 
minimize  quadratic  cost  due  to  error  between  aided  and  approximate  response. 
This  scheme  has  the  advantage  that  it  can  be  forced  over  the  bandwidth 
over  which  it  operates  the  closed-loop  controls. 

Single-Input  Single-Output  Transfer  Matching--Reference  40  presents  a 
method  which  has  been  shown  to  yield  excellent  results  on  particular  examples 
Very  favorable  comments  have  been  voiced  by  independent  users. 

Given  a  high  order  transfer  function 

1  +  a  s  + - +  a  s 

H(s)  =  K  - - - - -  (150) 

l+b,s+ - +  b  sn 

1  n 

where  m  <  n  which  is  to  be  approximated  by 
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L(s)  =  K 


1  +  C^S  + 


(151) 


c  s 
p 


1  +  d  s  t 

where 

p  <  q  <  n 

Define 

*<•>■*  iflg-l2 

The  X(u>)  is  expanded  in  a  MacLaurin  series.  This  leads  to  a  set  of  non¬ 
linear  algebraic  equations  for  determining  m  and  d^. 

Multiple-Input  Multiple-Output  Transfer  Matching — Reference  41  presents 
results  of  a  procedure  that  could  be  used  for  developing  low  order  matrix 
transfer  function  or  state  equation  approximations.  The  user  specifies 
(by  the  data  input)  the  bandpass  over  which  the  approximations  are  to  be  made. 

Reduced  Models  from  Noisy  Data — References  42  through  44  discuss  open- 
loop  low  order  approximations  in  a  noisy  environment.  They  approximate 
with  respect  to  noise  in  a  manner  analogous  to  frequency  domain  approxima¬ 
tions. 

Hierarchal  Control — Hierarchal  control  is  concerned  with  control  of  large 
systems.  The  flexible  aircraft  is  certainly  represented  as  a  large  system. 
Hierarchal  control  makes  synthesis  tractable,  in  part,  by  special 


: 


decomposition  techniques.  Recently,  two  treatises  have  been  published 
(References  48  and  49);  they  should  be  reviewed  for  applicability  to  the 
flexure  control  problem. 


TRUNCATION  OF  LARGE  MODELS  CONTAINING  RATES 


In  general,  the  physics  of  dynamical  processes  generate  a  model  in  the 
form  of 
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x2 
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(153) 


in  which  the  response  vector  contain  state  rates. 

In  this  case,  care  should  be  exercised  in  obtaining  the  truncated  model. 

The  truncated  model  can  be  obtained  either  before  or  after  differentiation 
with  different  results.  It  is  recommended  that  the  response  data  matrices 
C  ,  C  ,  C  ,  and  C  should  not  be  modified  before  the  truncation  operation. 

XX  1  ci  Z  l  zz 

Case  I:  Truncation  Before  Differentiation 


First  the  system  described  by  equations  (153)  is  truncated  to  obtain 


128 


and  then  substitution  for  x  is  made  to  obtain  the  truncated  state  model 


Case  II:  Truncation  After  Differentiation 


First  the  substitution  of  state  rates  into  response  equation  is  made 


and  then  Equation  (156)  is  truncated 


*i =  An  xi +  V 


r  -  (Cjj  +  C21  Ajj  +  C22  A21)  Xj 


(157) 


+  (C21  B1  +  C22  B2  +  D)  U 

As  can  be  seen,  the  truncated  models  represented  by  Equations  (155)  and 
(157)  differ  in  the  response  equations.  For  good  steady-state  response, 
the  differentiation  of  response  must  follow  the  reduction  process  as  presented 
in  Case  I.  In  the  residualization,  models  obtained  by  either  methods  are 
identical. 

In  Reference  4,  truncated  models  were  obtained  as  described  in  Case  I, 
before  differentiation  because  of  the  availability  of  C-5A  model  generating 
program,  generating  data  in  the  form  of  (153).  In  the  present  study  the 
response  data  in  the  form  of  (156)  was  available  only.  For  this  reason,  the 
responses  containing  state  rates  are  in  error  in  steady-state  for  truncated 
models. 

Effect  of  Truncation  on  Response  Rates 
Consider  the  system  described  by 
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11  12 


21  22 


—  — 

B, 

1 

1 

+ 

B0 

2 

2 

(158) 


r=[cn  Cl2] 


X. 


where  r  is  the  set  of  responses  whose  rates  are  required  as  design  responses. 
The  response  rates  can  be  obtained  either  before  or  after  truncation. 


Response  Rates  Before  Truncation- -Response  equation  is  differentiated  first. 

(159) 


ci2> 


Then  the  state  rates  are  substituted  from  (158)  and  the  resulting  equations 
are  truncated. 


This  yields 


*i =  An  xi +  Biu 


(160) 


r  =  (Cu  Au  +  C12  A21)  xx  +  <CU  Bj  +  CJ2  Bg)  u 
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Response  Rates  After  Truncation- -The  model  is  truncated  first  yielding 


*i =  An  xi +  V 


r  =  cu  xi 


(161) 


Then  the  truncated  response  equation  is  differentiated  and  the  truncated 
state  rates  are  substituted  from  (161)  yielding 


*■  =  CU41  ■  (C11  An>  xl  +  (C11B1)U 


(162) 


If  the  steady-state  value  of  the  response  rate  of  the  reduced  system  is 
required  to  be  zero,  then  the  truncated  model  represented  by  Equations 
(161)  and  (162)  should  be  used.  This  means  that  one  should  obtain  the 
response  rates  after  obtaining  the  reduced  model.  In  the  present  study, 
response  rates  were  obtained  before  truncation.  Because  of  this  some  of 
the  response  rates  in  the  truncated  model  do  not  satisfy  steady-state 
requirements. 

Wind  Simplification 

All  previous  discussion  has  been  concerned  with  simplifying  the  flexible 
body  representation.  This  is  the  primary  concern.  However,  time  and 
money  could  also  be  saved  by  using  the  simplest  wind  (and  wind  distribution) 
filter  that  met  requirements.  For  FLEXSTAB  C-5A  ALDCS  effort,  a  7th 
order  filter  was  used  for  both  synthesis  and  analysis.  It  is  believed  a  much 
lower  order  would  have  been  adequate  for  the  relatively  expensive  synthesis 
effort. 


ssBwrai 


In  summary,  a  cursory  survey  made  in  this  study  indicates  that  there  is  a 
dearth  of  closed-loop  and  a  surfeit  of  open-loop  methods  for  model  reduction. 
We  recommend  that  an  attempt  be  made  to  extend  the  closed-loop  scheme  pre¬ 
sented  in  Reference  32  into  an  engineering  tool  to  reduce  design  computation 
time.  Furthermore,  the  open-loop  methods  should  be  ranked  for  potential 
applicability  and  the  more  promising  refined  and  incorporated  into  KONPACT. 

C-5A  MODEL  COMPARISON  (FLEXSTAB  AND  HONEYWELL/GELAC) 

The  Honey well/GELAC  79th  order  model  (HG79)  for  the  cruise  flight  condition 
was  made  available  at  Honeywell  for  this  program.  The  42nd  order  model 
(HG42)  is  obtained  by  truncation  of  the  Wagner  dynamics.  Because  of  the 
way  the  unsteady  effects  were  modeled  in  the  79th  order  model,  the  42nd 
order  reduced  model  had  slight  errors  with  respect  to  the  original  Honeywell/ 
GELAC  42nd  order  model.  Table  15  shows  an  eigenvalue  comparison  between 
these  two  models  for  the  bending  modes.  Note  that  HG42  used  in  this  study 
has  less  damping  than  the  original  model.  This  required  more  iterations  on 
the  reduced  optimal  course  design. 

Table  15.  Root  Comparison  of  Free  A/C  with  HG42 
and  Original  HG42  Models 
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Air  Force  supplies  C-5A  data  for  the  cruise  flight  condition  in  the  form  of 
cards  (simulator  deck  data)  (  Reference  2,  5,  and  6).  This  data  was  con¬ 
verted  to  state  space  data  and  augmented  by  the  gust  model  and  actuators 
(used  in  the  Honeywell /Gelac  HG42  model).  The  resulting  system  was 
scaled  and  shuffled  to  bring  it  to  correspondence  with  Honeywell /GELAC 
data.  This  FLEXSTAB  model  (F42)  is  then  compared  and  correlated  to 
the  Honeywell /GELAC  model  (HG42). 

In  the  following  the  open-loop  comparisons  are  made. 

•  Steady-state  values 

•  rms  values  due  to  wind  gust 

•  Roots  (Eigenvalues) 

•  PSD  plots  due  to  wind  gust  and  white  noise 
inputs  to  aileron  and  inboard  elevator 

For  these  comparisons  the  FLEXSTAB  inboard  load  axis  system  (see  Figure  4) 
was  rotated  18.  6°  right  wing  tip  aft  to  be  comparable  to  the  HG42  model.  This 
rotation  was  performed  at  the  completion  of  this  study.  This  correction  is 
reflected  in  Tables  16,  17,  27b,  29  only  and  Figures  57,  59,  61,  63,  65,  and 
67  only. 

Table  16  compares  the  steady-state  1G  maneuver  responses  for  the  HG42 
and  F42  models.  The  rigid  body  responses  (w,  q,  a)  compare  very  favorably 
between  the  two  models.  Elevator  deflection  requirements  for  the  HG42 
model  are  50  percent  larger  than  those  for  the  F42  model;  this  is  to  be 
expected  since  the  HG42  model  data  have  been  corrected  for  the  disc  effects. 
The  bending  and  torsion  moment  data  also  agree  reasonably  closely.  The 
numerical  comparisons  between  the  bending  moments  and  torsion  moments 
are  very  good,  with  an  exception  being  that  of  T1  and  T2.  A  magnitude 
difference  is  to  be  expected. 
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Comparison  of  Free  A/C  with  HG42  and  F42 


The  torsional  moments  are  about  an  axis  parallel  to  the  wing  elastic  axis. 
The  magnitudes  of  the  torsional  moments  are  extremely  dependent  on  the 
chordwise  location  of  the  net  load  center  of  pressure  and  therefore  are 
extremely  math  model  dependent.  Sensors  outputs  (A21F,  A21R,  AFUS, 
and  TFUS)  are  exactly  the  same  as  should  be  expected.  The  first  hint 
that  there  are  major  differences  between  the  two  models  comes  in  viewing 
the  bending  mode  selection  data.  Magnitudes  and  signs  just  simply  do  not 
agree  any  place.  This  indicates  (but  does  not  prove)  that  there  may  be 
significant  difference  between  the  two  models. 

HG42  and  F42  model  free  aircraft  rms  responses  due  to  wind  gust  are 
compared  in  Table  17.  Generally,  the  comparisons  indicate  that  the  F42 
model  is  matching  the  HG42  model  well.  There  are  two  significant  dis¬ 
crepancies:  the  HG42  T2  torsion  moment  and  the  accelerometer  outputs 
A21F  and  A21R  are  larger  than  the  F42  data.  Since  exactly  the  same  gust 
field  is  "striking"  each  model,  data  indicate  that  the  HG42  model  is 
more  responsive  to  gust.  Other  than  that,  the  comparisons  are  very  good. 
This  would  indicate  that  the  F  model  is  a  good  approximation  to  the  HG 
model  for  the  wind  gust  inputs. 
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Table  18.  Root  Comparison  of 
F42  Models 
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In  Table  18  the  roots  of  HG79,  HG42  and  F42  aircraft  models  are  compared. 
The  root  values  are  quite  close.  They  are  good  enough  so  that  based  on 
these  data  alone,  one  would  expect  that  there  were  only  small  differences 
between  the  HG  and  F  model. 


Figures  56  through  7LTpresenl  the  power  density  plots  for  IIG42  model  and 

F42  model.  The  rms  values  as  obtained  from  PSD  plots  are  differ  ent  from 

« 

those  obtained  from  covariance  analysis.  This  difference  is  due  to  frequency 
factor  in  the  rms  values  computed  by  integrating  the  PSD  plots  as  given 
below. 
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The  PSD  plots  provide  two  valuable  services.  First,  they  are  an  excellent 
means  of  checking  models  to  ascertain  that  errors  have  not  been  caused 
by  the  data  handling  procedures.  Second,  the  qualitative  nature  of  the 
plots  provide  a  good  indication  of  where  the  problem  areas  are. 


The  PSD  plots  for  HG42  model,  as  presented  here,  compare  very  well 
with  the  previous  PSD  plots  obtained  during  previous  ALDCS  work 
(Reference  4).  This  shows  that  the  data  handling  in  KONPACT  with  the 
HG42  model  is  accurate. 


The  first  comparison  to  note  between  FLEXSTAB  and  Honeywell ,/GELAC 
models  is  to  compare  Figures  56  and  57  for  the  bending  moment  due  to 
wind  gusts.  For  the  HG42  model  resonant  peaks  occur  at  rigid  body 
near  0.2  cps  and  also  at  the  first  mode,  at  somewhat  less  than  1  cps. 


Figures  56  through  151  appear  behind  the  main  text  of  this  section 
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The  data  for  the  F42  model  show  a  lot  of  resonance  for  only  the  rigid  body. 
Similarly,  the  bending  moment  due  to  the  elevator  shows  for  the  HG42 
airplane  (in  Figuie  58)  has  resonant  peaks  at  the  rigid  body  first  and 
third  bending  modes.  The  F42  model  in  Figure  59  shows  a  majority  of 
the  resonance  at  the  rigid  body,  a  neglectable  amount  for  the  first  mode, 
and  then  apparently  picks  up  a  fourth  mode.  Figures  60  and  61  show  that 
the  HG42  model  has  significant  first  and  third  mode  energy;  the  F42  model 
has  significant  resonant  peaks  for  the  first  flexture  mode  only.  Torsional 
plots,  due  to  gusts  for  the  HG42  and  F42  models  in  Figures  62  and  63, 
show  that  the  HG42  aircraft  has  rigid  body,  first  and  third  flexure 
resonant  peaks;  the  F42  model  has  resonant  peaks  in  the  rigid  body  area 
only.  Torsion,  due  to  the  elevator,  is  shown  in  Figures  64  and  65.  For 
the  HG42  model  there  are  resonant  peaks  for  the  rigid  body,  the  first 
mode,  second  mode,  third  mode,  fourth,  sixth,  seventh  and  eighth  mode; 
for  the  F42  model,  there  are  resonant  peaks  for  the  rigid  body,  first 
mode,  fourth  mode,  and  seventh  mode;  these  latter  frequencies  are  very 
low  relative  to  the  values  that  we  see  on  the  HG42  model.  Torsion  due 
to  aileron  is  shown  in  Figures  66  and  67.  For  the  HG42  model,  there 
are  resonant  peaks  for  the  rigid  body,  first,  second,  third,  sixth,  seventh, 
and  thirteenth  mode;  for  the  F42  model,  there  is  also  an  additional 
resonant  peak  for  the  fifth  mode. 

The  third  flexure  mode  response  due  to  gusts  is  shown  in  Figures  68  and  69. 
The  1IG42  model  has  resonant  peaks  for  rigid  body,  first,  third,  and  fourth 
modes;  the  F42  model  has  resonant  peaks  for  the  rigid  body  only.  The 
third  flexure  mode  due  to  elevator  are  shown  in  Figures  70  and  71.  For  the 
HG42  model,  there  are  resonant  peaks  for  rigid  body,  third  and  fourth  modes 
and  the  F42  model  has  it  for  the  rigid  body,  first  and  third  modes.  For  the 


third  flexure  mode  bending  due  to  aileron,  plots  are  shown  in  Figures  72  and 
73.  There  is  qualitative  agreement  between  the  two  sets  of  data,  except 
that  the  rms  values  differ  by  a  factor  of  /2TT  .  Pitch  rate  responses  due  to 
gusts  are  shown  in  Figures  74  and  75.  The  comparisons  show  that  the  only 
responses  due  to  the  rigid  body  are  quite  good.  Pitch  rate  due  to  the  inboard 
elevator  are  shown  m  Figures  76  and  77.  The  comparisons  are  very  good. 
Pitch  rate  due  to  ailerons  are  shown  in  Figures  78  and  79  and  again  the 
comparisons  are  quite  good,  although  there  is  a  little  more  first  and  third 
mode  response  in  the  HG42  data. 

Further  open-loop  comparison  on  the  reduced  (residualized)  models  for 
FLEXSTAB  and  Honeywell/GELAC  is  given  in  the  next  subsection. 

RESIDUA  LIZA  TION  AND  TRUNCATION  STUDY  ON  FLEXSTAB  MODEL 

The  42nd  order  FLEXSTAB  model  is  used  to  obtain  three  different  models 
of  24th  order  for  residualization  and  truncation  study.  The  names  of  the 
models  and  the  procedure  to  obtain  them  are  shown  in  Table  19. 

Table  19.  Models  for  Residualization  and  Truncation  Study 
Obtained  from  the  F42  Model 


1  Model 

Procedure 

Residualize  states,  responses  and  sensors. 

F24RT 

Residualize  states  and  responses  and 
truncate  sensors. 

'  F24TT 

Truncate  states,  responses  and  sensors. 

Note  the  inboard  bending  moments  (Bl,  Tl)  (for  the  F24RR,  F24RT  and 
F24TT  models)  are  about  the  axis  system  shown  in  Figure  4. 


Tn  addition,  these  models  are  also  compared  with  HG24RR  obtained  from 
HG42  model.  First,  open-loop  (free  aircraft)  results  are  presented  and  then 
the  closed-loop  (ALDCS)  results  are  presented. 

The  following  performance  measures  are  used  for  comparison: 

«  Steady -state  values 
»  rms  values  due  to  wind  and  pilot 
*  Roots  (Eigenvalues) 

c  Time  response  plots. 

Open-] mop  Results 

Table  20  compares  the  steady-state  free  aircraft  1G  maneuver  responses 
for  the  F24RR,  F24RT,  and  F24TT  models.  The  quantities  w,  a,  q, 
flexure  modes,  and  bending  moments  compare  favorably  among  all  models. 
Table  21  illustrates  the  effect  of  residuaiization  and  truncation  on  the  steady- 
state  value  of  sensors.  (Note  that  truncation  here  means  truncation  after 
differentiation.)  In  Reference  4,  there  was  much  closer  correspondence 
among  the  RR,  RT,  and  TT  models.  (There,  truncation  was  done  before 
differentiation. )  Also  note  that  the  steady-state  value  of  T5  for  the  F24TT 
model  is  about  20  times  lower  than  that  for  the  other  models.  (This  again 
is  due  to  truncation  after  differentiation  operation. )  This  fact  was  not  known 
before  this  study  was  undertaken. 

One  other  significant  thing  to  note  is  that  the  steady -state  values  of  the  bending 
and  torsion  moment  rates  for  F24TT  model  are  not  zero  as  should  be  expected. 
For  this  reason  the  respor.se  rates  should  be  obtained  after  the  reduced 
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Table  20.  Steady-State  Comparison  of  Free  A/C  with  F24RR,  F24RT 
and  F24TT  Models  for  1G  Maneuver  (Continued) 


Table  20.  Steady-State  Comparison  of  Free  A/C  with  F24RR,  F24RT 
and  F24TT  Models  for  1G  Maneuver  (Concluded) 
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model  is  obtained.  (In  this  study  the  response  rates  were  obtained  for  the 
F42  model  and  then  the  reduced  models  were  obtained. ) 

Table  21.  Free  A/C  1G  Responses 


'V"''\vModels 

Sensors'''''.. 

F42 

F24RR 

F24RT 

F24TT 

0.04377 

0.04377 

0.04378 

0.04326 

AFUS,  g's 

-1.00000 

-1.00000 

-0.78000 

-0.75000 

A21R,  g's 

-1.00000 

-1.00000 

-3.90000 

-3. 60000 

Table  22  shows  the  comparison  of  free  aircraft  rms  responses  due  to  gust 
input  for  HG24RR,  F24RR,  F24RT,  and  F24TT  models.  As  indicated 
by  the  rms  comparisons  between  HG42  and  F42  models,  the  rms  comparisons 
between  HG24RR  and  F24RR  are  also  very  good.  This  means  that  the  different 
reduction  procedures  do  not  affect  the  rms  values  significantly  (as  they  do  for 
the  steady-state  values). 

Similar  comparisons  are  made  for  the  pilot  response  for  the  HG24RR, 

F24RR,  F24RT,  and  F24TT  models  as  shown  in  Table  23,  Again  the  com¬ 
parisons  among  the  models  would  indicate  that  the  FLEXSTAB  model  is  good 
for  design  of  controllers  to  take  care  of  the  rms  disturbances  (i.  e,,  to 
handle  the  control  problem,  at  least  at  low  frequencies). 


In  Table  24,  the  roots  of  HG24RR,  F24RR,  and  F24TT  models  for  free  A/C 
are  compared.  Again,  the  root  values  are  close.  The  damping  ratios  of 
the  F24RR  model  (as  shown  in  Table  24  and  in  Table  18  for  the  F42  model) 


Table  24.  Root  Comparison  of  Free  A/C  with  HG24RR,  F24RR 
and  F24TT  Models 


HG24RR 

F24RR 

F24TT 

Association 

Real 

Real 

Real 

U> 

n 

c 

“n 

C 

e 

0.  a 

1.  6663 

0.4981 

1. 5501 

0.5785 

1.  6987 

0.5854 

V  *1 

5.9309 

i - 

0.0854 

5.  5283 

0.0892 

5.5501 

0.0921 

V  ^  2 


V  ^  3 


V  % 


V  ^  5 


V  \ 


T\r 


eo 


A21RL 


GLAF 


MLC1 


Pilot  Filter 


Gust  Filter 


11.4585  0.0207 


14.5012  0.0286 


16.0682  0.0447 


17.4781  0.0235 


19.1533  0.0252 


-6.0000 


-7. 5000 


-7. 5000 


-6.  0000 


-0.0200 


-0.0100 


-0. 1000 


-0. 2100 


15.  5911  0.0699 


17.2662  0.0491 


18.3701  0.  0245 


19.3531  0.0315 


22.  1171  0.0487 


-6.0000 


-7. 5000 


-7.5000 


-6.0000 


-0.0200 


-0.0100 


-0. 1000 


-0. 2100 


15.5416  0.0707 


17.4777  0.0531 


18.3842  0.0263 


19.3590  0.0323 


22.0976  0.0499 


-6.0000 


-7. 5000 


-7. 5000 


-6.0000 


-0.0200 


-0.0100 


-0. 1000 


-0. 2100 


Table  24.  Root  Comparison  of  Free  A/C  with  HG24RR,  F24RR 
and  F24TT  Models  (Concluded) 


HG24RR 


F24RR 


F24TT 


Association 

Real 

Real 

Real 

a) 

n 

C 

% 

c 

0) 

n 

Gust  Filter 

-0.2100 

-0. 2100 

-0. 2100 

Wing  Kussner 

-9. 1560 

-9. 1560 

-9. 1560 

1st  Order 
Delay 

-13.4270 

-13. 4270 

-13. 4270 

2nd  Order 
Delay 

9. 8020 

0. 8165 

9. 8020 

0. 8165 

9. 8020 

Tail  Kussner 

-18.  4930 

-18.4930 

-18.4930 

change  very  little  due  to  model  reduction  process.  Previous  experience  had 
indicated  that  there  would  be  typically  about  50  percent  more  damping  ratio 
on  the  flexture  modes  for  the  lower  order  residualized  models  than  for  the 
high  order  models. 

Figures  80  througn  99  contain  the  time  response  plots  for  HG24RR  model  and 

F24RR  model  to  elevator  input.  The  elevator  input  to  produce  this  series  of 

plots  is  shown  in  Figures  96,  97,  98,  99;  this  is  the  amount  of  elevator 

deflection  required  to  pull  a  1-g  incremental  maneuver.  The  responses 

(Bl,  Tl,  a,  q,  7]  ,  7|  7]  fi  ,  6  .,  *  )  for  the  two  models  are  all  very 

1  d  6  ei  ei  eo  J 

similar. 


Closed-Loop  Results 

Table  25  compares  the  steady-state  ALDCS  aircraft  1G  maneuver  responses 
for  the  F24RR,  F24RT,  and  F24TT  models.  (Note:  ALDCS  controller 
is  designed  with  F24RR  model  as  explained  in  Section  V  and  the  same 
controller  is  used  with  the  models  F24RT  and  F24TT).  The  following 
(Table  26)  illustrates  the  comparison  of  steady-state  values  of  sensors. 

The  different  model  reduction  procedures  seem  to  have  less  effect  in 
closed-loop  performance  measures. 

By  comparing  the  data  for  the  free  aircraft  from  Table  20  and  the  data 
from  ALDCS  aircraft  in  Table  25  we  get  the  data  as  shown  in  Table  27a 
and  27b. 

The  maneuver  relief  for  the  F24RR  model  for  B1  is  nearly  the  exact  value 
specified  (namely  0.  7).  The  F24RT  model  gets  a  little  more  bending  relief 
than  required  and  the  F24TT  model  shows  that  only  about  4/5  required 
bending  relief  is  achieved.  Bending  reliefs  achieved  with  different  models 
is  comparable  to  the  data  that  were  obtained  during  the  ALDCS  design 
phase  using  the  Honeywell/GELAC  model,  (Reference  4). 

The  second  line  presents  comparable  data  for  torsion.  Here  it  is  seen  that 
the  relative  torsion  values  for  the  F24RR,  F24RT,  and  F24TT  models  are 
0.2983,  0.3116,  and  0.0761.  These  lorsional  moment  ratios  were  produced 
relative  to  the  axis  system  shown  in  Figure  4.  If  the  SS  loads  of  Tables 
20  and  25  are  rotated  18.  6°  to  agree  with  the  Honeywell/GELAC  model, 
then  Table  27b  results.  These  data  agree  very  favorably  with  the  Honeywell/ 
GELAC  in  that  the  torsional  moments  increase  with  the  ALDCS  aircraft. 
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Table  25.  Steady-State  Comparison  of  ALDCS  A/C  with  F24RR,  F24RT  and 
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Table  25.  Steady-State  Comparison  of  ALDCS  A/C  with  F24RR,  F24RT  and 
F24TT  Models  for  1G  Maneuver  (Concluded) 
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Table  26.  Closed-Loop  1G  Responses 


^\Models 

Sensors"''. 

F24RR 

F24RT 

F24TT 

TFUS  rad/seo 

0.04375 

0.04824 

0.04916 

AFUS.g's 

-1.00100 

-0. 89800 

-0. 89400 

A21R,  g’s 

-1.00100 

-1. 12800 

-1. 11700 

Table  27a.  Steady-State  Maneuver  Load  Relief  Ratio 
(Moment  ALDCS/Moment  Free) 

(Figure  4  Load  Axis  System) 


F24RR 

F24RT 

F24TT 

B1 

0.  6991 

0.7635 

0.  7821 

T1 

0.2983 

0.3116 

0.0761 

Table  27b.  Steady-State  Maneuver  Load  Relief 
(Moment  ALDCS/Moment  Free) 
(Figure  3  Load  Axis  System) 


F24RR 

F24RT 

F24TT  j 

B1 

0.6686 

0.7289 

0.7328 

T1 

1.8046 

2.009 

2.1702 

In  Table  28  steady-states  are  compared  for  the  free  and  ALDCS  aircraft 
using  the  HG24RR  and  F24RR  models.  Table  29  compares  bending  and 
torsion  data  relief  for  the  F24RR  and  HG24RR  models.  It  is  seen  that 


eady-State  Comparison  of  Free  A/C  and  ALDCS  A/C  with  F24RR 


'KM>  tfi  fcPS  k.]h  U>  *  *4  .tjui  a  *  u\  o| 

♦  *  V  *  *  ♦  ♦  ?iT  ?  ”  M  n**  *"ncl'i 

1  b!  l-l  HI  U  < ■  •  U.I  U]  ijjlli  ill  i,*l  UJ  l! J L* I  J,  J,  u  J.l 


°  -*  If*  e-'L»  K  *4  m  m  ,o  u.  IV.  n  oln^  *4  &  K.  rr  e>  □  u  m  □  f 

r\j  k>  n  n.  ojk.  w  w  u'  ai  .t  it.  ir  tr  coin  k  k  *  «-*  <\j  o  in  L  o>  t 

i2  *c  *  *2  rfe ,ft  n  +  *•  +  °  K  ^  ‘0p. «,  ^  £  *,  ~  S  £  E  S I 


i  i  i 

UJ  UJ  lu 

7>  *0  W 

c*  O  a 

Nlft  (M 


nriooc  n»imr  i  nunc  o  o  u 

+  **t**»*±i  +  *^*  +  l  |  l 

UJ  U'  (•'  b'  U  I|J  'l'  I  .1  I  J  I  111  til  b‘  >'  (l  tu  UJ  Ul 

W  4  WK  .1  OM  H  t’  r  O'  r .  •-«  IT  (Ti  IT-  K. 

a'U'cr^-o'  re  m  a.-  m  a  m  ^  n  m  s  r  <y  o> 

u«i«wr  l/i  O'  rl  ri  I  in  IA  i£  s.  u>  in  H  N 


in  -  in  ,7  -<  m  . 

f  O  •<  L3  n  H  O  1 

•  »•••!> 

IU  lu  U  tv  In  lu  Itl  l 
■A  M  o  n.  ji  in  r,  * 
»-i  r  &  t*  j  >q  < 

c  to  j  o  <u  rn  n 


ft/rf  niT  ro  ^  k.  »j.i  u>m  m  -Hui  aj  ^  ~t  r>  -*U.  in  w 


DN-!®r*n'«N^®MNN  D  I/I  l/l  *f>  0  IT  .  n  ^ 

>ooo  o-oe^ooxoo  a  o  o  w  9  o  o  f-  :<  o  » 

>«»*«|l**|4*||««*  |«  .11  I  I  I  I  I  I 

iiu  w  u  ui  y  w  uj  uj  uj  w  uj  uj  uj  jujw  ujij  uj  uj  uj  *.  uj  u 

O  kll/l‘4MTM\l"4C  nc«c  CT  CT  o  7  .»  U 

n-  *i  <y  —  ^-n 

►  t>nrvjcorvor~  —  ciyr~t\j^  -  <7  l/i  —  7  ^  >c  -  r- 

DN4  4nMlMNM(V,ni|A>h4A<  <C  V  n  a.  a.  -  a,- 


®  *  n*  r*  —  ir*  u>  e.  .7  it  —  7  ir  c  >r  ir 

o  ©  —  o  «  o  ©  o  —  o  o  —  a  o  —  •:>  o 

. . . 

UJ  UI  UJ  UJ  U  Uj  Ui  Uj  UI  U  UJ  Uj  U:  UJ  U  U.  UJ 

<r  r  4  tr  airin  p*  -  % 

’j  0*  j  .t  •>  iO  4 
..  l\J'C'4-r-*j-.1tT4n7«-r^ 
—  a  —  j-f^cro'— -i/*.  or  —  c*  a  j  ^ 


mo  on  n  n  |o  o  mu 

UJ  uj  UJ  uj  UI  uj  (uj  ui  uj  m 

UI  l/l  1/1  1/1  1/1  1/1  Il/l  l/l  VM/> 

VN.  VN  s\  K\  s\ 


In  icn  p  to 
p  ar  j  _j  rp  jr 

fl  •*  »  •  4|l  | 


U  U  f  U  M 

cd  a  lu  tr  is  in  m  m  u.'  o  ra  u,  o  tct 

!'-J_ll'i->-JI/--J_ll/i_J.Jl/>_J  — J* 
1  I  •  V  I  I  \  I  ■  \  I  I  s  l  ih 


.  TPT  w  T  r  M  r  T’k-*r-r*~.-r*r-»-f-'»-’Y--rT--r-rT--r-r-*'f'*- 

lOlinrinn  i)  fninT  L'nn  HMiuiurMomi  nno 

- 7rtT7*liT7Z7r7T?7TrTTTT 

..(I'M  MU'  M  MA'MI-<IVMI-«l»4i*^l-«k-il-ll-lh-ih-IMW<t~lMMfc-l 
_ _ _ ^ 


tf  OUA  f 
M  M  U  f 
ft'UH’V 

o  *-  -  o  * 


’  7  Z7  7.  T  X  at 

I  Ul  III  I<|  it  In  111  UI 

■  r  sr  r  r  t  r  r 

i  o  or*  on  no 

ryrrurrvtruu 


■»  •*  Z  m  n  : 
.  7  r*  r<  7  r 
M  M  O  ►-«  ► 

'  f  j  mj  o  i 


o  f-  »v 
J77U 
Z  U  UJ  IS 

z  r 

UJ  O  o  ill 

ur  r*. 


►-  *-  *H  K  *-  Z  *-  I 

7  Z  «f  Z  2  H  7  . 

U  UI  UJ  UI  U  l  I 

Z  X  t  7  r  ! 

o-o  no  o  i 

a.  r  r  z  r  : 


!wi*7wi*7  m  m 
inrrr 70mm 
40MI-40MW00 
:  _i  n  r/i  »  o  (/i  _j  _j 
‘  hi  7  V  1.17  IV  hi  hi 

>>u»o>ujO>> 

»  fit  ►*  r 


1“  »1  Z  <V 

W  Z  «? 

H  H  t 
IPMM 
J2D!W 

n  m  n  « 

U.  jf  rfl. 


(!)  7  J  U  7  J  II  : 
7  0Ui?Pli!7  ( 
ynOMHOHi 
O  CM  O  CM  O  l 

z  ^  z  a?  z  - 

III  O  III  O  lit  I 

if,  u  «»  ui  w.  «*  r  i 


i  o  r»  <  r 

AVI*  V 
:  OQIUP 


UWUt-l*.  UUJWUiftU- 

rin  c.  o  «i  d  c  u  cnti«i 

c  iu  r  o  m  e  r  ui  n  c<  ui  u.1 


v  (**  iv  v  (v'l/v  /y  jy  r 


Hiv  H'iviv  /v[<v  y  (V  »v  ivl/v  (V  v  r.  ' 


the  torsion  moment  for  the  HG  model  has  increased  as  was  mentioned  before. 
However,  the  magnitude  of  this  increase  is  greater  than  that  obtained  during 
the  past  ALDCS  design  study.  This  difference  might  be  due  to  the  fact 
that  the  ALDCS  controller  designed  is  not  a  very  refined  one.  It  is  seen 
that  the  T). *s  are  markedly  different  between  the  HG24RR  and  F24RR  models 

Table  29.  Steady-State  MLC  for  HG24RR  and  F24RR  Models 
(Figure  3  Load  Axis  System) 


F24RR 

HG24RR 

El 

0.6686 

0.7002 

T1 

1.8046 

1.7428 

_ 

and  this  raises  a  question  on  the  accuracy  of  the  F  model  at  higher  frequencies. 
(Th>  discrepancy  was  also  observed  while  comparing  the  HG42  and  P’42  models. ) 

Table  30  shows  the  comparison  of  ALDCS  aircraft  rms  responses  due  to 
gust  input  for  IIG24RR,  F24RR,  F24RT  and  F24TT  models.  Table  31 
shows  the  relative  rms  relief  as  obtained  from  Tables  22  and  30. 

An  observation  is  that  the  different  improvement  obtains u  between  the  RR, 

Rl,  TT  models  would  lead  one  to  markedly  expect  the  different  results  de¬ 
pending  upon  which  of  the  models  were  used.  Experience  in  the  past  during 
the  design  of  the  ALDCS  was  that  the  A  LDCS  rms  relief  in  gust  was  the  same 
regardless  of  the  model  used.  This  difference  is  because  of  the  different 
truncation  procedures  used  in  this  study. 


Table  30.  RMS  Comparison  of  ALDCS  A/C  vdth  HG24RR,  F24RT  and 
F24TT  Models  for  1  FPS  Disturbance  on  T) 


A  second  observation  is  in  regard  to  the  torsion  relief.  During  the  ALDCS 
design  it  was  observed  that  torsion  increases  could  be  maintained  below 
5  percent.  The  torsion  reductions  shown  here  for  the  FLEXSTAB  models 
are  again  attributed  to  the  rotated  axis  system  (cf.  Figures  3  and  4). 


Table  31.  Relative  ALDCS  RMS  Relief  in  Gusts 


HG24RR 

F24RR 

F24RT 

F24TT 

B1 

0. 5032 

0.6045 

0.  8757 

0.7775 

T1 

1.3071 

0.4820 

0.6528 

1.3451 

Table  32  presents  the  rms  comparisons  for  the  ALDCS  aircraft  due  to 
pilot  disturbance,  and  Table  33  presents  the  results  of  a  combined  dis¬ 
turbance.  Table  34  compares  the  rms  responses  for  the  free  and  ALDCS 
A/C  with  the  F24RR  model. 

Roots  for  the  ALDCS  aircraft  for  all  the  four  models  (HG24RR,  F24RR, 
F24RT,  and  F24TT)  are  compared  in  Table  35.  Here  the  very  low  damp¬ 
ing  rate  ratio  is  generated  on  the  third  flexure  mode  for  the  HG24RR 
model.  To  improve  the  damping  ratio  it  was  necessary  to  put  low  pass 
filter  into  the  elevator  channel.  As  will  subsequently  be  shown,  the  elevator 
drives  the  third  mode  rather  hard  on  the  HG  model.  The  table  does 
show  that  the  rigid  body  objectives  for  all  four  aircraft  are  as  desired. 

Table  36  compares  the  roots  of  the  free  and  ALDCS  aircraft  with  the  F24RR 
model.  The  ALDCS  aircraft  displays  better  root  positions  than  do  2S  the 
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Table  32.  RMS  Comparison  of  ALDCS  A/C  with  HG24RR,  F24RR,  F24RT  and 
F24TT  Models  for  1  FPS  Disturbance  on  1] 
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Table  34.  RMS  Comparison  of  Free  A/C  and  ALDCS  A/C  with  F24RR  Model  for 
1  FPS  Disturbance  on  1)  and  T| 


iSftglfoilr 


free  aircraft.  Note  that  the  damping  ratio  for  the  F24RR,  F24RT,  and  F24TT 
models  in  the  third  bending  mode  is  increased.  This  is  not  what  would  have 
been  expected  from  previous  experience  in  designing  the  HG  model. 

The  time  response  plots  for  the  closed-loop  systems  are  presented  in  Figures 
100  through  151.  The  following  closed-loop  models  are  considered. 

•  SAS  aircraft  HG24RR 

•  SAS  aircraft  F24RR 

•  ALDCS  aircraft  HG24RR 

•  ALDCS  aircraft  F24RR 

•  ALDCS  aircraft  F24RT 

•  ALDCS  aircraft  F24TT. 

It  is  shown  in  Figures  100  and  101  that  the  SAS  improves  the  pitch  rate 
response  for  both  the  HG24RR  and  F24RR  models.  The  basic  C-5A  with 
SA.S  is  considered  to  have  very  desirable  handling  qualities  and  ALDCS  is 
designed  to  provide  the  same  qualities  for  pitch  rate. 

Figures  102,  103,  104,  and  105  show  that  the  pitch  rate  response  for  the 
ALDCS  aircraft  is  very  close  to  the  SAS  aircraft  previously  displayed  and 
show  that  the  control  design  technique  met  their  objectives.  For  the  response 
shown  for  the  ALDCS  HG24RR  aircraft  (in  Figure  102),  it  is  apparent 

that  there  is  a  considerable  amount  of  the  third  flexure  mode  appearing  in  the 
pitch  rate  response;  during  the  previous  ALDCS  design  (Reference  4), 
a  better  looking  pitch  rate  response  was  obtained  for  the  ALDCS  aircraft. 

This  mildly  oscillatory  portion  of  the  responses  is  due  to  the  low  damping 


Table  36.  Root  Comparison  of  Free  A/C  and  ALDCS 
A/C  with  F24RR  Model 


F  roe 

ALDCS 

Heal 

Real 

Association 

U) 

n 

C 

0) 

n 

C 

9.  a 

1. 5501 

0. 5785 

1. 7393 

0.7354 

V  T 

1*  l 

5. 5283 

0. 0892 

3. 7410 

0.2501 

T2*  f2 

3*  3 

15. 5911 

0.0699 

14. 2402 

0. 1852 

17. 2662 

0.0491 

17. 2567 

0.0169 

T)  -n 

5'  5 

18. 3701 

0.0245 

18. 3744 

0.0216 

• 

Ti6'  t,6 

19. 3531 

0.0315 

19. 1853 

0.0324 

V  f7 

22. 1171 

0.0487 

21. 9848 

0.0569 

«a 

-6.0 

-1. 2552 

®ei 

-7.  5 

-8.0200 

1 

sco 

-7.5 

-7.5 

A21RL 

-6.0 

30. 5312 

0.8669 

GLAF 

-0.  02 

30. 5312 

0.8669 

MLC1 

-0.01 

-0.01 

Pilot  Filter 

-0.  1 

-0.  1 

Gust  Filter 

-0.  21 

-0.  21 

Gust  Filter 

-0.  21 

-0.  21 

Wing  Kussner 

-9. 156 

-9. 156 

1st  Order  Delay 

-13.427 

-13.427 

2nd  Order  Delay 

9.  802 

0. 8165 

9.  802 

0.  8165 

Tail  Kussner 

-18.493 

-18. 493 

For  the  bending  moment  response,  there  are  two  objectives.  The  first  is 
that  for  the  steady-state  value  of  the  bending  moment  with  the  ALDCS 
aircraft  should  be  70  percent  of  that  over  the  free  aircraft.  Second,  it  is 
very  desirable  that  there  be  as  little  overshoot  as  possible  in  the  bending 
moment  response  because  each  overshoot  represents  needless  damage 
to  the  aircraft.  The  ALDCS  designs  under  this  contract,  j  shown  in 
Figures  106,  107,  108,  and  109,  all  display  rather  good  bending  moment 
response;  there  is  little  overshoot.  The  major  objection  to  the  response 
is  that  on  the  HG24  model  (Figure  106)  there  is  the  third  mode  oscillation 
which  would  be  eliminated  in  a  more  refined  design;  for  example,  the  third 
bending  mode  filters  could  be  put  into  attenuate  the  elevator  effects  on  the 
third  mode,  HG24RR  responses  exhibit  a  little  more  of  the  third  mode  than 
would  be  desired. 

The  torsion  moment  responses  for  the  ALDCS  aircraft  are  detailed  in 
Figures  110,  111,  112,  and  113.  As  with  the  bending  moment,  the  torsion 
responses  should  be  smooth  with  just  a  little  overshoot.  The  results  for 
the  HG24RR  model  are  quite  good  except  for  the  third  mode;  this 
effect  could  be  eliminated  in  a  more  refined  design  procedure  than 
was  accomplished  here.  The  results  for  the  F24RR  model  again  look 
surprisingly  close  to  those  for  the  HG  model  with  a  lot  of  the  third  bending 
mode  appearing.  The  results  for  F24RT  and  F24TT  models  are  very 
different  from  what  we  expected  to  see.  The  reason  could  be  the  truncation 
after  differentiation  procedure  that  was  adopted  during  the  study.  The 
results  for  the  torsional  version  show  that  the  general  objectives  for  the 
ALDCS  are  achieved. 
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First  flexure  mode  responses  for  the  ALDCS  aircraft  are  shown  in  Figures 
114,  115,  116,  and  117,  Note  that  the  first  flexure  mode  for  the  HG24RR 
model  has  ve*y  low  damping.  During  the  previous  ALDCS  design  a  higher 
damping  wa«  achieved.  The  F24RR  model  ALDCS  results  are  qualitatively 
satisfactory  although  the  deflections  are  large  relative  to  that  found  on  the 
HG24  model.  The  results  for  the  F24RT  and  F24TT  models  appear  again 
to  be  too  different  than  previous  experience  would  indicate. 

The  transient  responses  for  the  Liird  bending  mode  for  the  ALDCS  aircraft 
are  shown  in  Figures  118,  119,  120,  and  121.  For  the  HG24RR  ALDCS 
response,  the  third  bending  mode  has  little  damping;  this  checks  with  the 
damping  ratio  of  0.0051  that  was  shown  previously.  This  corresponds  to 
an  amplitude  ratio  of  0.966  per  cycle.  With  the  HG34  model  as  indicated 
previously,  the  third  mode  is  very  difficult  to  control.  It  takes  more 
effort  and  more  tailoring  to  suppress  this  mode.  The  other  ALDCS  traces 
also  show  a  well  damped  third  mode  oscillation. 

The  sixth  flexure  mode  results  for  the  ALDv_o  aircraft  are  shown  in  Figures 
122,  123,  124,  and  125.  The  sixth  flexure  mode  shows  rather  modest 
damping.  In  viewing  these  results,  it  should  be  noted  that  the  amplitudes 
and  the  responses  are  quite  small;  the  normalized  plots  that  we  are 
presenting  here  tend  to  over-emphasize  the  importance  of  these  terms. 

Elevator  response  plots  for  the  ALDCS  aircraft  in  Figures  126  through  137 
show  considerable  activity  which  reflects  the  feedback  from  the  flexure 
modes  into  the  ALDCS  controller.  The  aileron  response  plots  are  shown 
in  Figures  138  through  145.  For  ALDCS,  the  aileron  response  should  be 
compared  with  the  bending  moment  response.  The  aileron,  to  a  high  degree 
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establishes  bending  relief.  Therefore,  the  aileron  response  needs  to  be 
phased  in  such  a  manner  that  it  prevents  the  overshoot  in  bending  moment 
response.  The  ALDCS  objective  of  getting  the  bending  moment  response 
to  look  like  a  second  order  system  with  a  damping  ratio  of  about  0.  7  can  be 
achieved  by  phasing  the  aileron  correctly  with  the  feedback  system.  For 
example,  for  the  HG24RR  results,  it  is  seen  that  the  bending  moment  has 
a  flat  spot  between  1  and  2  seconds  and  then  peaks  at  around  2.  7  seconds 
and  then  backs  off  by  about  5  percent;  these  are  results  in  Figure  106. 

The  aileron  transient  response  is  shown  in  Figure  138.  It  peaks  at  about 
2.0  seconds  and  stays  there  until  about  0.2  seconds  and  retrenches  at  about 
5  percent.  If  the  aileron  response  was  reshaped  so  that  the  portion  between 
0.7  of  a  second  and  1.  6  seconds  had  the  outer  end  moved  from  1. 6  to  about 
2.  7,  the  bending  moment  response  of  Figure  106  would  have  been  improved 
because  of  the  elimination  of  the  overshoot.  In  any  event  since  the  step 
command  response  is  really  a  worse  case  situation,  the  response  for  normal 
input  with  the  system  as  shown  here  would  have  very  little  overshoot  and 
would  be  very  desirable.  Comparisons  of  the  ALDCS  results  w'ill  shov; 
that  the  aileron  phasing  is  quite  close  to  that  desired  in  order  to  prevent  the 
bending  mode  excesses.  (Note  in  Figures  126  to  130,  BDEI  is  equal 
to  DELEI)  . 

The  angle  of  attack  response  plots  are  shown  in  Figures  146  through  151. 

In  summary,  it  can  be  said  that  on  the  basis  of  these  comparisons,  the 
FLEXSTAB  produces  valid  vehicle  models  for  preliminary  design  of  CCV 
control  laws. 
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Figure  56.  PSD  Plot  of  Free  A/C  with  HG42  Model  (Bl/ETAG) 
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Figure  57.  PSD  Plot  of  Free  A/C  with  F42  Model  (Bl/ETAG) 
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Figure  58.  PSD  Plot  of  Free  A/C  with  HG42  Model  (Bl/UDELEI) 
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Figure  61.  PSD  Plot  of  Free  A/C  with  F42  Model  (Bl/UDELA) 
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Figure  62.  PSD  Plot  of  Free  A/C  with  HG42  Model  (Tl/ETAG) 
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Figure  66.  PSD  Plot  of  Free  A/C  with  KG42  Model  (Tl/UDELA) 


Figure  68.  PSD  Plot  of  Free  A/C  with  HG42  Model  (UE3/ETAG) 
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Figure  (0.  PSD  Plot  of  Free  A/C  with  HG42  Model  (UE3/UDELEI) 
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Figure  72.  PSD  Plot  of  Free  A/C  with  I1G42  Model  (UE3/UDELA) 
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Figure  74.  PSD  Plot  of  Free  A/C  with  HG42  Model  (Q/ETAG) 
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Figure  81.  Command  Response  Plot  of  Free  A/C  with  F24RR  Model  (Bl) 
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Figure  82.  Command  Response  Plot  of  Free  A/C  with  HG24RR  Model  (Tl) 
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Plot  of  Free  A/C  with  HG24RR  Model  (ALPHA) 


2  to* 


6tO* 


I 

soo  •«* 

X 

Sifi  *t 

I 

S 26  •£ 

I 

S88  4  £ 

■ 

$98  T 

I 

S08  *£ 

X 

592  *  £ 

X 

522  ’£ 

I 

$89  *  £ 

I 

$99  4£ 

I 

509  't 

I 

$95  '£ 

I 

$25  T 

T 

$«*  *X 

5 

$99  4t 

i 

$09  4£ 

i 

$9C4f 

i 

S2r  4£ 

i 

S82  T 

i 

STT 

i 

S02  T 

i 

$91  4  £ 

i 

S2T  4 £ 

i 

S80  4£ 

i 

S*»0  4E 

i 

SCOT 

i 

S96  42 

i 

526  *2 

i 

SOO  *2 

i 

ST  *2 

i 

$C8  *2 

i 

592  *2 

i 

$22  *2 

i 

$89*2 

i 

599  4  2 

i 

SC9  4 2 

i 

£95  4  2 

i 

525*2 

i 

$P9  4 2 

i 

ST  *2 

i 

SD9  *2 

i 

59£  *2 

i 

$2£  *2 

i 

£82  *2 

i 

ST  *2 

i 

Stz  4 2 

i 

£91  *2 

i 

S2T  *2 

I*-—'* 

i 

$80  *2 

oo 

i 

590  4Z 

Q 

i 

soo  *z 

z: 

i 

$96*1 

o 

i 

$26  4 1 

C_) 

i 

588  4 1 

LU 

i 

$**B4I 

oo 

i 

SC9  M 

- - 

i 

592  4  T 

i 

S22  T 

Ixl 

i 

S89  4 T 

s; 

i 

ST  4  T 

i— t 

i 

$09  4  T 

t— 

i 

$  95  4  I 

i 

$  25  4  T 

i 

SS*>  •  T 

i 

£♦»•»•  T 

i 

$  09  4T 

i 

$9f  T 

i 

SZE  4 T 

i 

5  82  4 T 

i 

5*2  4  I 

i 

♦  32  4T 

i 

S9T  4 T 

i 

S2T  4 T 

i 

$80  4T 

i 

$90  4t 

i 

SCO  4 T 

i 

$96  4 

i 

SZ6  4 

r 

$88  4 

i 

$98  4 

i 

£06  4 

i 

$92* 

i 

S?2  4 

i 

$89  4 

i 

$99  4 

i 

SG9  4 

i 

$95  4 

i 

$25  4 

X 

589  4 

I 

ST  4 

I 

S09  4 

I 

S9C  4 

I 

S2C  4 

I 

SBZ  4 

•  I 

S92  4 

•I 

S02  4 

•I 

S9I  4 

*1 

S2t 4 

•1 

$*fl  * 

«x 

S9I  4 

•  X 

SOT 

— X 

(PBU) 


203 


,  ^  1..  I--.-'.  ■  ■  .  , 


Figure  85.  Command  Response  Plot  of  Free  A/C  with  F24RR  Model  (ALPHA) 
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Figure  86.  Command  Response  Plot  of  Free  A/C  with  HG24RR  Model  (Q) 


Figure  87.  Command  Response  Plot  of  Free  A/C  with  F24RR  Model  (Q) 
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Figure  88.  Command  Response  Plot  of  Free  A/C  with  HG24RR  Model  (ETA1) 
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Figure  89.  Command  Response  Plot  of  Free  A/C  with  F24RR  Model  (ETA1) 
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Figure  90.  Command  Response  Plot  of  Free  / 
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Figure  92.  Command  Response  Plot  of  Free  A/C  with  HG24RR  Model  (ETA6) 


Figure  93.  Command  Response  Plot  of  Free  A/C  with  F24RR  Model  (ETA6) 
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Figure  96.  Command  Response  Plot  of  Free  A/C  with  HG24RR  Model  (DELEI) 
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Figure  97.  Command  Response  Plot  of  Free  A/C  with  F24RR  Model  (DELEI) 
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Figure  99.  Command  Response  Plot  of  Free  A./C  with  F24RR  Model  (DELEO) 
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Command  Response  Plot  of  ALDCS  A/C  with  HG24RR  Model  (Q) 
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Figure  104.  Command  Response  Plot  of  ALDCS  A/C  with  F24RT  Model  (Q) 
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Figure  106.  Command  Response  Plot  of  ALDCS  A/C  with  HG24RR  Model  (Bl) 
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Figure  107.  Command  Response  Plot  of  ALDCS  A/C  with  F24RR  Model  (Bl) 
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Figure  109.  Command  Response  Plot  of  ALDCS  A/C  with  F24TT  Model  (Bl) 
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Figure  110.  Command  Response  Plot  of  ALDCS  A/C  with  HG24RR  Model  (Tl) 
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Figure  111.  Command  Response  Plot  of  ALDCS  A/C  with  F24RR  Model  (Tl) 
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Figure  112.  Command  Response  Plot  of  ALDCS  A/C  with  F24RT  Model  (Tl) 
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Figure  113.  Command  Response  Plot  of  ALDCS  A/C  with  F24TT  Model  (Tl) 
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Command  Response  Plot  of  ALDCS  A/C  with  F24RT  Model  (ETA1) 
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Figure  119.  Command  Response  Plot  of  ALDCS  A/C  with  F24RR  Model  (ETA3) 
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Figure  120.  Command  Response  Plot  of  ALDCS  A/C  with  F24RT  Model  (ETA3) 
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Figure  124.  Command  Response  Plot  of  ALDCS  A/C  with  F24RT  Model  (ETA6) 
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Figure  125.  Command  Response  Plot  of  ALDCS  A/C  with  F24TT  Model  (ETA6) 


Figure  126.  Command  Response  Plot  of  ALDCS  A/C  with  HG24RR  Model  (DELEI) 
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Figure  130.  Command  Response  Plot  of  ALDCS  A/C  with  HG24RR  Model  (DELEIDOT) 
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Figure  131.  Command  Response  Plot  of  ALDCS  A/C  with  F24RR  Model  (DELEIDOT) 
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Figure  132.  Command  Response  Plot  of  ALDCS  A/C  with  F24RT  Model  (DELEIDOT) 
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Figure  133.  Command  Response  Plot  of  ALDCS  A/C  with  F24TT  Model  (DELEIDOT) 


Figure  135.  Command  Response  Plot  of  ALDCS  A/C  with  F24RR  Model  (DELEO) 
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Command  Response  Plot  of  A LDCS  A/C  with  F24RT  Model  (DELEO) 
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Figure  137.  Command  Response  Plot  of  ALDCS  A/C  with  F24TT  Model  (DELEO) 
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Figure  139.  Command  Response  Plot  of  ALDCS  A/C  with  F24RR  Model  (DELA) 
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Figure  142.  Command  Response  Plot  of  ALDCS  A/C  with  HG24RR  Model  (DELADOT) 
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Figure  144.  Command  Response  Plot  of  ALDCS  A/C  with  F24RT  Model  (DELADOT) 


Figure  146,  Command  Response  Plot  of  SAS  A/C  with  HG24RR  Model  (ALPHA) 
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Figure  147.  Command  Response  Plot  of  SAS  A/C  with  F24RR  Model  (ALPHA) 


Figure  148.  Command  Response  Plot  of  ALDCS  A/C  with  HG24RR  Model  (ALPHA) 
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SECTION  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  objectives  of  this  study  were  threefold:  1)  develop  an  interface 
program  between  FLEXSTAB/LSA  and  Air  Force-owned  optimal  control 
programs  (DLAK  and  FFOC),  2)  demonstrate  the  interface  by  using  C-5A 
model  and  Active  Lift  Distribution  Control  System  Design  Procedure,  and 
3)  conduct  a  brief  study  on  model  reduction  procedures  for  design.  The 
major  emphasis  has  been  on  the  interface  software  development. 

These  objectives  were  primarily  met.  The  algorithms  used  in  the  interface 
and  the  results  of  demonstration  example  are  documented  in  this  report. 

The  developed  programs  are  listed  in  AFFDL-TR-75-146,  Volume  II. 

Users  information  for  the  programs  is  given  in  AFFDL-TR-75-146  Volume  III. 

In  the  following,  the  results  and  recommendations  for  future  studies  pertaining 
to  the  work  in  the  area  of  analysis  and  synthesis  and  software  developments 
are  presented. 

SIGNIFICANT  RESULTS 

•  The  work  reported  here  established  the  total  software  system 
approach  to  Active  Control  Technology  (ACT)  and  Control 
Configured  Vehicles  Design  (CCV)  problems. 
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e  The  chief  benefit  of  the  program  was  to  provide  integrated  software 
for  rapid  development  of  vehicle  model  and  control  laws  for  ACT 
and  CCV  studies. 

•  The  results  of  demonstration  example  with  C-5A  ALDCS  show  that 
FLEXSTAB  program  system  generates  valid  models  for  preliminary 
control  system  synthesis, 

•  The  results  also  show  that  the  residualization  procedure  is  better 
than  truncation  procedure  when  higher  order  models  are  reduced 
to  low  order  models  for  design. 

RECOMMENDATIONS  FOR  FUTURE  ANALYTICAL  WORK 

•  Development  of  efficient  model  reduction  procedures  is  needed 
for  feedback  design  of  high  order  models  representing  unsteady 
aerodynamics  and  bending  modes. 

•  Development  of  design  procedures  using  residualized  Riccati 
equations  via  singular  perturbation  and  boundary  layer  techniques 
is  needed  to  reduce  design  cost  for  high  order  systems. 

•  Development  of  efficient  solution  procedure  is  needed  for  constrained- 
optimal  control  problems  to  enforce  conventional  design  criteria 

into  optimal  active  control  synthesis  techniques. 

RECOMMENDATIONS  FOP.  FUTURE  SOFTWARE  DEVELOPMENT  WORK 


•  Development  of  geometry  definition  interface  module  is  needed 


to  generate  vehicle  configuration /performance  sensitivities  for 
Control  Configured  Vehicles  design. 

•  Automatic  weight  selection  and  response  selection  modules  should 
be  developed  to  speed  up  the  quadratic  design  process. 

•  Algorithms  developed  for  the  Constrained-Optimal  Control  problem 
should  be  implemented. 

•  Existing  optimal  control  routines  should  be  revised  foi  higher 
computational  efficiency. 

•  Optimal  multirate  digital  control  analysis  and  design  capabilities 
should  be  added  using  the  Air  Force-owned  programs. 

CONCLUSIONS 

A  large-scale  software  system  for  automatic  modeling  and  optimal  design 
of  active  control  systems  was  developed  in  this  study.  The  software  system 
was  used  to  design  controller  for  the  C-5A  vehicle  with  Active  Lift  Distribu¬ 
tion  Control  System  design  procedure.  Automatic  dimensioning  and  user 
oriented  statement  features  considerably  improved  the  interface  data  mechanics 
for  optimal  design. 


REFERENCES 


1.  A.  J.  VanDierendonck  and  G.  L.  Hartmann,  "Quadratic  Methodology,  " 
Volume  II  -  Documentation  for  Computer  Program,  Honeywell  Report 
No.  F0161-FR,  October  1973, 

2.  G.  R.  Hank  el  al. ,  "A  Method  for  Predicting  the  Stability  Character¬ 
istics  of  Control  Configured  Vehicles,  "  Volume  II  -  FLEXSTAB 
2.01.00  User's  Manual,  APFDL-TR-74-91 ,  November  1974. 

3.  A.  F.  Konar  and  J.  K.  Mahesh,  "Computer  Programs  for  Active 
Control  Technology,  "  Volume  II,  KONPACT  Program  Listings, 
AFFDL-TR-75-146  Vol,  I,  June  1976. 

4.  R.  B.  Beale,  M.  F.  Barrett,  and  C.  R.  Stone,  "Control  Design  for 
the  C-5A  Active  Lift  Distribution  Control  System,"  Honeywell  Report 
No.  RF-09523C,  Volumes  I  and  n,  April  1974. 

5.  R,  Schwanz,  C.  Stockdale,  "Effect  of  Structural  Idealization  on 

the  Aeroelastic  Stability  and  Control  Parameters  of  the  C-5A  Aircraft,  " 
AFFDL-TM-75-122-FGC,  December  1975. 

6.  C.  Stockdale  and  G.  Grimes,  "A  Formulation  of  Dynamic  Load  Equations 
using  the  Level  2  FLEXSTAB  Data  Base  for  Control  Configured  Vehicles 
Design,  "  AFFDL-TM-75- 120,  February  1976. 

7.  Phillips,  A.  C.,  Rosen,  R.  W. ,  "A  Note  on  Dynamic  Data  Storage  in 
Fortran  IV,  "  Computer  Journal,  Vol.  18,  No.  4,  November  1975. 


8.  M.  F.  Barrett  and  R.  C.  McLane,  "Design  Procedure  and  Results 

for  ALDCS,  "  MR  12223,  S&RC  Honeywell,  Minneapolis,  September  1973. 

9.  E.  N.  Tinoco  and  J.  E.  Mercer,  "FLEXSTAB,  A  Summary  of  the 
Functions  and  Capabilities  of  the  NASA  Flexible  Airplane  Analysis 
Computer  System,  "  NASA  CR-2564,  December  1975. 


273 


REFERENCES  (Continued) 


10.  R.  D.  Miller  et  al. ,  "Feasibility  of  Implementing  Unsteady  Aero¬ 
dynamics  into  the  FLEXSTAB  Computer  Program  System,  "  NASA 
CR  132530,  October  1974. 

11.  A.  J.  VanDierendonck  et  al. ,  "Application  of  Practical  Optimal 
Control  Theory  to  the  C-5A  Load  Improvement  Control  Systems 
(LICS),"  AFFDL-TR-73-122,  October  1973. 

12.  A.  F.  Konar,  "Development  of  Weapon  Delivery  Models  and  Analysis 
Programs,  "  AFFDL-TR-71-123,  Volume  I,  April  1972. 

13.  G.  L.  Hartmann  et  al. ,  "F-8C  Digital  CCV  Flight  Control  Laws,  " 

NASA  CR  2629,  February  1976. 

14.  M.  S.  Borow  et  al. ,  "Navy  Digital  Flight  Control  System  Development,  " 
Honeywell  Document  No.  21857-FR,  GAP,  Minneapolis,  December 
1972. 

15.  C.  R.  Abrams,  "A  Performance  Index  for  Response  Evaluation  of 
Highly  Augmented  Military  Aircraft,  "  NADC-AM-7103,  Naval  Air 
Development  Center,  Warminster,  Pa.,  12  October  1971. 

16.  "Military  Specification  -  Flying  Qualities  of  Piloted  Airplanes,-" 
MIL-F-8785B  (ASG),  7  August  1969. 

17.  "Military  Specification  -  Control  and  Stabilization  Systems:  Automatic, 
Piloted  Aircraft,  General  Specification  For,  "  MIL-C-18244A  (WEP), 

1  December  1962. 


18.  L.  Edinger  and  T.  Lahn,  "C-5A  Data  Base  for  Load  Alleviation  and 
Mode  Stabilization  Program,  "  Report  20564-DB1,  Honeywell  Inc., 
Government  and  Aeronautical  Products  Division,  Minneapolis,  Minn., 
1  April  1968. 


....  Ji 


274 


REFERENCES  (Continued) 


"Aircraft  Load  Alleviation  and  Mode  Stabilization  (LAMS):  C-5A 
System  Analysis  and  Synthesis,  "  Technical  Report  AFFDL-TR-68-162, 
Wright- Patterson  Air  Force  Base,  Ohio,  November  1969. 

C.  R.  Stone,  M.  D.  Ward,  C.  A.  Harvey,  M.  E.  Ebsen,  E.  E.  McBride, 
and  W.  W.  Hollenbeck,  "Studies  on  the  Compatibility  of  Relaxed 
Static  Stability  and  Maneuver  Load  Control  to  C-5A-Type  Aircraft,  " 
Volume  I,  Technical  Report  AFFDL-TR-72-38,  Wright- 
Patterson  Air  Force  Base,  Ohio,  June  1972. 

Edinger,  Schenk,  and  Curtis,  "Study  of  Load  Alleviation  and  Mode 
Suppression  (LAMS)  on  the  YF-12A  Airplane,  "  NASA  CR-2158  , 

June  1972. 

"Electronic  Load  Improvement  Control  System  (LICS)  Cost  and 
Effectiveness  Study  for  the  C-5A,  A  Report  to  the  C-5A  Independent 
Structural  Review  Team,"  13  September  1972. 

G.  Stein,  A.  H.  Henke,  "A  Design  Procedure  and  Handling-Qualities 
Criteria  for  Lateral- Directional  Flight  Control  Systems,  "  AFFDL- 
TR-70-152,  May  1971. 

A,  J.  VanDierendonck,  "Design  Method  for  Fully  Augmented  Systems 
for  Variable  Flight  Conditions,"  AFFDL-TR-71-152,  January  1972. 

A.  J.  VanDierendonck,  "  Duadrater  Methodology :  A  Snort  Course 
on  the  Application  of  Quadratic  Optimal  Control  Theory  to  the  Design 
of  Practical  Flight  Control  Systems  (or  any  other  similarly  Modeled 
System),  "  prepared  for  AFFDL  under  contract  No.  F33615-72-C-2008 . 

C.  R.  Stone,  M.  D.  Ward,  C.  A.  Harvey,  M.  E.  Ebsen,  E.  E.  McBride, 
and  W.  W.  Hollenbeck,  "Studies  in  the  Compatibility  of  Relaxed  Static 
Stability  and  Maneuver  Load  Control  to  C-5A  Type  Aircraft,  " 
AFFDL-TR-72-38,  Volume  II,  June  1972. 


275 


REFERENCES  (Continued) 


27.  C.  R.  Stone,  "From  LICS  to  ALDCS  with  F.  C.  37  of  the  C-5A,  " 
Honeywell  Research  Memo,  MR  12184,  2!)  June  1973. 

28.  C.  R.  Stone,  "LICS  to  ALDCS  Revised,"  Honeywell  Research  Memo, 

MR  12195,  10  July  1973. 

29.  C.  R.  Stone,  "Modeling  of  Flexible  C-5A's  for  ALDCS,"  Honeywell 
Customer  Engineering  Letter  to  Lockheed  Georgia  Company,  C.  E.  L. 
No.  ALDCS-2,  27  July  1973. 

30.  J.  G.  Truxal,  "Automatic  Feedback  Control  System  Synthesis," 
McGraw-Hill  Book  Co.  Inc.,  New  York  1955. 

31.  A.  F.  Konar  and  J.  K.  Mahesh,  "Digital  Control  Systems  for  Tactical 
Fighters,  "  AFFDL-TR-73-119,  Volume  I,  June  1974. 

32.  Petar  V.  Kototovic  and  Ricard  A.  Yackei,  "Singular  Perturbation  of 
Linear  Regulators:  Basic  Theorems,"  IEEE  Trans.  Auto.  Ctrl,, 
Volume  AC-17,  No.  1,  February  1972,  pp.  29-37. 

33.  Richard  A.  Yackei  and  Petar  V.  Kokotovic,  "A  Boundary  Layer 
Method  for  the  Matrix  Riccati  Equation,  "  IEEE  Trans.  Auto.  Ctrl. , 
Volume  AC-18,  No.  1,  February  1973,  pp.  17-24. 

34.  Y.  Shamask,  "Stable  Reduced-Order  Models  Using  Pade-Type 
Approximations,  "  IEEE  Trans.  Auto  Ctrl.,  October  1974,  pp.  615-616. 


35.  Y.  Shamask,  "Linear  System  Reduction  Using  Fade  Approximation  to 
Allow  Retention  of  Dominant  Modes,  "  Int.  J,  Ctrl. ,  Volume  21,  No.  2, 
1975,  pp.  257-272. 


36.  Y.  Shamask,  "Model  Reduction  Using  the  Routh  Stability  Criterion  and 
the  Fade  Approximation  Technique,  "  Int.  J.  Ctrl. ,  Vol.  21,  No.  3, 
1975,  pp.  475-484. 


_ - 


276 


REFERENCES  (Continued) 


37.  M.  Lai  and  R.  Mitra,  "A  Comparison  of  Transfer  Function  Simplifica¬ 
tion  Methods,  "^EEEjrrarisi_Auto:_CtrL_,  October  1974,  pp.  617-618. 

38.  Maurice  F.  Hutton  and  Bernard  Friedland,  "Routh  Approximations 
for  Reducing  Order  of  Linear,  Time -Invariant  Systems"  IEEE  Trans. 

Auto.  Ctrl. ,  Vol.  AC-20,  No.  3,  June  1975,  pp.  329-337. 

39.  E.  E.  Yore  and  Y.  Takahashi,  "identification  of  Dynamic  Systems  by 
Digital  Computer  Modeling  in  State  Space,  "  Journal  of  Basic  Engineering, 
June  1967,  pp.  295-299. 

40.  T.  C.  Hsia,  "On  the  Simplification  of  Linear  Systems,  "  IEEE  Trans. 

Auto  Ctrl.  June  1972,  pp.  372-374. 

41.  R.  3.  Beale  and  N.  E.  Miller,  "Turbine  Engine  Control  Synthesis, 

Volume  III:  Experimental  Engine  Identification  and  Modeling,  "  Aero 
Propulsion  Lab  Technical  Report  AFAPL-TR-75-14,  March  1975. 

42.  R.  Genesio  and  P.  Pome,  "Identification  of  Reduced  Models  from  Noisy 
Data,  "  Int.  J.  Ctrol,  Volume  21,  No.  2,  1975,  pp.  203-211. 

43.  T.  Subba  Rao,  "An  Innovation  Approach  to  the  Reduction  of  the  Dimensions 
in  a  Multivariable  Stochastic  System,  "  Int.  J„  Ctrl.,  Volume  21,  No.  4., 
1975,  pp.  673-680. 

44.  Lewis  Meier  and  David  G.  Luenberger,  "Approximation  of  Linear 
Constant  Systems,"  IEEE  Trans.  Auto.  C'rl.,  October  1967,  pp.  585- 
588. 

45.  C.  T,  Leondes  and  John  F.  Yocum,  "Optimal  Observers  for  Continuous 
Time  Linear  Stochastic  Systems,  "  Automatica,  Volume  II,  1975, 

pp.  61-73. 


277 


REFERENCES  (Concluded) 


46.  Robert  G.  Schwendler  and  Richard  H.  MacNeal,  "Optimum  Structural 
Representation  in  Aeroelastic  Analyses,  "  Flight  Dynamics  Laboratory 
Report  ASD-TR-61-680,  March  1962. 

47.  Mihajlo  D.  Mesarovic,  Theory  of  Hierarchal  Multilevel  Systems, 

Academic  Press,  1970. 

48.  Krzysztof  Cichocki,  "Multilevel  Control  Systems,  "  University  of 
Minnesota,  May  1975  (primarily  a  translation  of  Reference  49). 

49.  W.  Findersen,  "Multilevel  Control  Systems,  "  (in  Polish),  PWN, 

Warzawa,  1974. 

50.  "A  Method  for  Predicting  the  Stability  Characteristics  of  Control 
Configured  Vehicles,  "  Volume  I  -  FLEXSTAB  2.01.00  Theoretical 
Description,  Boeing  Commercial  Airplane  Company,  AFFDL-TR-74-91, 

November  1974. 

/ 

51.  A.  F.  Konar  and  J.  K.  Mahesh,  "Computer  Programs  for  Active 
Control  Technology,  "Volume  III,  User's  Manual,  AFFDL-TR-75-146, 

Vol.  IB,  June  1976. 

52.  "New  Short  Period  Handling  Quality  for  Fighter  Aircraft,  "  Boeing 
Document  D6-17841  T/N,  October  1965. 

53.  Kisslinger,  R.  L. ,  and  Wendl,  M, ,  "Survivable  Flight  Control  System 
Interim  Report  No.  1,  Studies  Analysis  and  Approach,"  AFFDL-TR-71-20 
Supplement  1,  May  1971. 

54.  Dornfeld,  G. ,  Fraser,  R.  ,  "Review  of  AFFDL  Program,  AFLOADS," 

The  Boeing  Commercial  Airplane  Company,  Contract  F33615-75-C- 
3132  Management  Review  Report,  August  1976. 


278 

•fr  U.  S.  GOVERNMENT  PRINTING  OFFICE:  1976  —  657-630/1057 


turn 


■■  ■  .Mw-dto* 


